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Abstract. We recall the form factors f^ N corresponding to the A-extension C(N,N; X) 
of the two-point diagonal correlation function of the Ising model on the square lattice and 
their associated linear differential equations which exhibit both a "Russian-doll" nesting, and 
a decomposition of the linear differential operators as a direct sum of operators (equivalent 
to symmetric powers of the differential operator of the complete elliptic integral E). The 
scaling limit of these differential operators breaks the direct sum structure but not the 
"Russian doll" structure, the "scaled" linear differential operators being no longer Fuchsian. 
We then introduce some multiple integrals of the Ising class expected to have the same 
singularities as the singularities of the n-particle contributions x to the susceptibility of 
the square lattice Ising model. We find the Fuchsian linear differential equations satisfied 
by these multiple integrals for n = 1,2,3,4 and, only modulo a prime, for n — 5 and 6, 
thus providing a large set of (possible) new singularities of the x ■ We get the location 
of these singularities by solving the Landau conditions. We discuss the mathematical, as 
well as physical, interpretation of these new singularities. Among the singularities found, 
we underline the fact that the quadratic polynomial condition 1 + 3w + Aw 2 = 0, that 
occurs in the linear differential equation of actually corresponds to the occurrence 

of complex multiplication for elliptic curves. The interpretation of complex multiplication 
for elliptic curves as complex fixed points of generators of the exact renormalization group 
is sketched. The other singularities occurring in our multiple integrals are not related to 
complex multiplication situations, suggesting a geometric interpretation in terms of more 
general (motivic) mathematical structures beyond the theory of elliptic curves. The scaling 
limit of the (lattice off-critical) structures as a confluent limit of regular singularities is 
discussed in the conclusion. 
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1 Introduction 

This paper displays a selection of works and results that have been obtained by the authors 
in collaboration with B.M. McCoy, W. Orrick and J. -A. Weil. It also provides new ideas and 
viewpoints at the end of Subsection 13.41 in Section [5] and in the Conclusion [TU1 We also give 
new results of linear differential operators modulo prime that had not been published before in 
Appendix IC.3I and Appendix IC. 41 

The two dimensional Ising model in zero magnetic field is, historically, the most important 
solvable model in all of theoretical physics. The free energy [91], the partition function on 
the finite lattice [67] and the spontaneous magnetization [921 1125| were computed long ago by 
Onsager, Kaufman and Yang. These computations, and subsequent studies of the correlation 
functions |68[ 1123] , form the basis of scaling theory and of the renormalization group approach 
to critical phenomena. 

Let us first recall the form factors [25] of the lattice Ising model. Our starting point will be 
the expansions of the diagonal correlations in an exponential form |123j . both for T < T c 

C.(N,N) = (1 -t)V*exp (E F S) M 

with 

t = (smh(2E v /k B T)smh(2E h /k B T)^j 2 
and for T > T c 

oo / oo \ 

C + (N, N) = (l- t)V* £ G^+D exp £ FiS >JV+1 (1.2) 

n=0 \n=l / 

with 

t = ((sinh(2E v /k B T) sinh(2E h / k B T)Y 

where E h and E° are the horizontal and vertical interaction energies of the Ising model. We 
will restrict in the following to the isotropic Ising model. For diagonal correlation functions, 
there is no difference between the isotropic and anisotropic models: the diagonal correlations 
are functions of the modulus k = smh(2E v /k B T) smh.(2 E h /k B T). The difference comes with 
off-diagonal modes and is sketched in [26j . 

When the exponentials in (II. ip and (|1.2j) are expanded, the correlations can also be written 
in what is called a "form factor" expansion 

CL (N, N) = (l- t) 1 / 4 (l + fj , (1.3) 

oo 

C + (N,N) = (l-t)^^f% n N +1 \ (1.4) 

n=0 



The form factor fffl N is interpreted as the "j-particle" contribution to the two-point correlation 
function. It is natural to consider A-extensions |83t 1123] of the previous functions 

C- (N, N; A) = (1 - t)V4 h + fjA 2 ™/^ , (1.5) 
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C + (N, N ; A) = (1 - t)V* A 2n /& +1) (1.6) 

n=0 

(7) 

which weight each fj^ N by some power of A, and to interpret A as being analogous to a coupling 
constant in a quantum field theory expansion. Such A-extensions naturally emerge from the 
Fredholm determinant framework in [123] . We will present new integral representations for F^^j, 

G^m^ and fff N in Section [2j We will see that they are much simpler, and more transparent, 
than the forms obtained from C(M,N) of |123| by specializing to M = N. 

The diagonal correlations C(N, N) have the property, discovered by Jimbo and Miwa [60] in 
1980, that their log-derivatives are solutions of the "sigma" forrrfl of a Painleve VI equation 

. _.d 2 a\ 2 „ T , /, ^.da \ 2 da f , .da 1 \ / da \ „, 

where cr is defined for T < T c as 

, 1N din C-(N,N) t 
a N (t) = t(t - l) 

with the normalization condition 

C-(N,N) = l + 0(t) for t->0 (1.8) 
and, for T > T r , as 



dlnC + (JV,iV) 



with the normalization condition 

C + (N,N) = ^^t N / 2 (l + 0(t)) for t-0, (1.9) 

where (a) at = T(a + N)/T(a) denotes the Pochhammer symbol. 

One can easily verify that (|1.7p . the A^-dependent sigma form of Painleve VI, is actually 

covariant by the Kramer s-Wannier duality 

(t,a,a',a") (j,j,a-ta',t 3 a" 

On another hand, Jimbo and Miwa introduced in |60j an isomonodromic A-extension of 
C(N,N). Remarkably this more general function C(N,N;X) also satisfies \27\ 193] the Pain- 
leve VI equation (|1.7p . The motivation of introducing an isomonodromic parameter A, in the 
framework of isomonodromy deformations, is, at first sight, quite different from the "coupling 
constant Fredholm-expansion" motivation at the origin of the form factor A-extensions (11.51) 
and (jl.6p . In [27] we have shown that these two A-extensions are actually the same by demon- 
strating that the recursive solutions of (jl.7p . analyticll in t 1 / 2 , agree with (11. 5j) and (jl.6p where 
the /w jy's are obtained from C±(N, N; A), the A-extension of C±(N, N). The normalization 
condition (jl.8p fixes one integration constant in the solution to (jl.7p . We find that the second 
integration constant is a free parameter, and, denoting that parameter by A, we find that our 
one-parameter family of solutions for C_ (JV, iV) can be written in a form structurally similar 
to the right hand side of (11. 5p . Furthermore, we have confirmed, by comparison with series 



(i) 

expansions of the multiple integral formulas for fw N derived in Section [21 that this family of 
solutions is, in fact, identical to C—(N,N;X) as defined in (11. 5ft . Similarly, the condition (11. 9p 
gives rise to a one-parameter family of solutions for C+ (N, N) that is identical to (II. 6|) . 



1 We use a variable t which is the inverse of the one of Jimbo and Miwa [60] . 

2 The A-extensions (|1.5p and (|1.6p are analytic at t ~ in t for T <T C and, when T > T c , analytic in t for iV 
even, and in f 1 / 2 for odd. 
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2 New integral representations for the /jy ]v' s 

The form factor expressions for the two-point correlation functions C(M, N) of |87l [88l [Ml [93"1 
112311124] are obtained by expanding the exponentials in (jl.ip . and (jl.2p . in the form given in [123] 
as multiple integrals, and integrating over half the variables. The form of the result depends 
on whether the even, or odd, variables of |123j are integrated out. For the general anisotropic 
lattice, one form of this result is given, for arbitrary M and N, in [93]. When specialized to the 
isotropic case the result is 

where s denotes sinh(2-fiT), and where 



I] h ik) (H x n) cosMV^J (2.1) 

l<i<k<j I \n=l / \ n=l / 



with 



s + ^ — cos 4> n — ^ + ^ — COS (j) n ) - 1 

2 \ V2 



1/2 



1 ,y . , 2(x^ fc ) 1/2 sin((^-^)/2) 

smh7 n = s H cos n - 1 , h ik - 



T iij I I 7 — tn, 

S J J 1 X{X k 

For T < T c , let us first recall equation (3.15) of Wu's paper [122] . which reduces, for the 
diagonal correlations C(N,N), to 

-1/2 



x / d?C N ((1 " aaO(l - « 2 /0) V2 77^2 > (2-2) 



(1 - r 1/4 C(iV, N )~ l + J^yJ d & N ((! " «20(1 " «2/e)) 

ce'-o 2 ' 

where 02 is i 1//2 . Comparing with (jl.3p we see that the second term in (|2.2p is j^n = ^NN- 

Performing the change of variables £ = z\ and £' = 1 /^2 ; deforming the contour of integration 
for both z\ and Z2 (one has to consider only the discontinuity across the branch culd running 
from to 02), and rescaling z\ and Z2, in, respectively, x\ = zi/a.2 and X2 = Z2/CH2, we obtain: 



n Jo Jo 

(1 - tXiX 2 ) 



x 2 (l - Xi){l - txi)J 

Similarly, when T > T c , the leading term for G\^ N is given by equation (2.29) of [122] 



'N,N 

1 



p(i) _ r a) _ - 1 / dz 

W " ^ " 2- y o ((1 _ ,1/2 Z)(1 _ ,1/2,-1)) V2 



3 For T < T c , a 2 = t 1/2 < 1 
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which, after deforming the contour of integration to the branch cut, and scaling z = t^^x, 
becomes 

n Jo 

= M^ 2Fl ^N + \-N + l-ty (2.3) 

where 2 Fi(a, b; c; z) is the hypergeometric function [57] . 

When the low temperature expansion of Section 3 of Wu [122] is performed to all orders, we 
find that ([TT]) holds with 



F 



(2n) 



1 ffr zfd Zj " / (l-a 2 z 2j )(l-a 2 /z 2j ) V /2 



N,N 



TT *j TT _ 

n (2vr)« 7 11 1 - 2^ i+a 11 V (1 - a 2 z 2j -i)(l - "2/^-1) 



from which, after deformation of integration contours and rescaling, one obtains, for T < T c , 
the following new integral representation of Fj^y(t): 



*NM T ) nn 2n llL 1 _ tx . x . +1 Ll\ X2 . (1 _ X2 . l)(1 _ tx2 . l) • W 



{ _ 1)n+H n(N + l) ,1 2n n . _ ^ _ ^ y/2 



J'=l 



Similarly for T > T c the expansion of Section 2 of Wu [122] is performed to all orders and 
we find that CC2D holds with given by and 



n+l ! 2n 

(2n+l) 



r 1 

N,N 



n+l n 

x n ((1 - a 2 - 1 %i-i)(l - a, 1 /^-!))" 172 II (C 1 " « 2 " 1 ^)(1 - «2 V^)) 1/2 • 
i=i i=i 

Changing variables and deforming contours, we obtain: 

£7V(2n+l)/2+2n <•! 2ra+l -, 2n 



n+l , v 1/2 n 

n( (.-^-^-,) ) n«i-^(i-^/^. 



X 



The form factor expressions are then obtained by expanding the exponentials. Thus we find, 
for T < T c , that the form factors in (jl.5p read 

(niy TT 2n J ax xx ^.(1 - X 2j _l)(l - tX 2 j-{)) 

n n 

X I 1(! - tx 2j -ix 2h )~ 2 (a?2j-i - Z2fc-i) 2 (£2j - x 2 k) 2 , (2.5) 

j=l fc=l l<j<k<n 

and, for T > T r , the odd form factors in (11.61) read 



f (2n+l),.x _ ,((2n+l)AT/2+n(n+l)) _ 1 
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x 



r i 2n+l n+1 

J Y[ x k dxk Yl ((! ~ x 2j)(l - tx 2 j)x 2 j) 1/2 



k=i j=i 

n+1 n+1 n 

X ((1 - X 2 j-l)(l ~ tX2j_l)x 2 j_l)~ 1/2 J ^(1 - tX 2 j-iX 2k )~ 2 

j=l j=l k=l 

x (a;2i-i - x 2 k-if (x 2 j-x 2k ) 2 , (2.6) 

l<i<fc<n+l l<j<k<n 

where the last product in (|2,6p has to be taken to be equal to unity for n = 0, 1. We note that the 
factors l/(n!) 2 and l/(n!(n + 1)!) in (|23J) and QZJft, arise because the integrands are symmetric 
functions of the variables x 2 j and x 2 j-±, separately. This is to be contrasted with (|2.ip . where 
there is no separation in the odd and even integrals . In the simplest case the previous integral 
representation (|2.6p gives f^ N (t) defined by where one recognizes the Euler representation 
of an hypergeometric function. 

Do note that the (G^^, N+l ) decomposition in (jl.2p is not unique. In contrast, the 

form factor expressions (12.5p . (|2.6j) are unique and well-defined. 

3 Fuchsian linear differential equations for fppjy{t) 

(i) 

We use formal computer algebra to study the functions f^ N . We obtain the Fuchsian linear 

(i) 

differential equations satisfied by the fj^ N for fixed j < 9 and arbitrary N. We also find the 

truly remarkable result that the f^N^ and Inn are ea ch solutions of linear differential oper- 
ators which have a nested "Russian-doll" structure. Beyond this "Russian doll" structure, each 
linear differential operator is the direct sum of linear differential operators equivalent^ to sym- 
metric powers of the second order differential operator corresponding to fjP N , (or equivalently 
to the second order differential operator Lg, corresponding to the complete elliptic integral E). 
A direct consequence is that the form factors f^^ l \ and /jy jy- are polynomials in the complete 
elliptic integrals of the first and second kinds, K and E 

K = 2 iq(l/2,l/2;l;t), E = 2 iq (1/2, -1/2; 1; t) . (3.1) 

(2) 

A simple example is /q = K(K — E)/2. 

In previous studies on the Ising susceptibility [126|. 1127} 11281 1129] , efficient programs were 
developed which, starting from long series expansions of a holonomic function, produce the 
linear ordinary differential equation (in this case Fuchsian) satisfied by the function. In order 
for these programs to be used to study the /^^y's we need to efficiently produce long (up to 
several thousand terms) series expansions in t of the /j^^y's. We have done this by use of both the 

integral representations (|2.5p . (|2.6p and the representations of /jy jy in terms of theta functions 
of the nome of elliptic functions, presented in [93] . 

We have obtained the Fuchsian linear differential equations satisfied by the (diagonal) form 

(i) 

factors /-ft jy up to j = 9. The analysis of these linear differential operators shows a remarkable 
Russian-doll structure similar to the nesting of (the differential operators of) the x s found 
in [1261 1127} 1128} I129j . Specifically we find that the expressions f]p N , /jyjy) /i/jyi frpN are 
actually solutions of the linear ODE for /jy jy, and that fffl N , /jyjyj /jy jy> /jy jy are actually 

(8) ' ' . ' 

solutions of the linear ODE for f^ N- ^ n addition, we find that all the linear differential operators 



4 For the notion of equivalence of linear differential operators, see |108l 11141 1115] . 
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(i) 

for the fftftS have a direct sum decomposition in operators equivalent to symmetric powers of 

the linear differential operator corresponding to fff N - Consequently, all the /jyjy's can also be 
written as polynomials in terms of the complete elliptic integrals E and K. The remainder of 
this section is devoted to the presentation of these results. 

3.1 Fuchsian linear differential equations for f^^iy 1 ^ 

The linear differential operator Fg(N) which annihilates f^ N has the following factorized form 

F 9 (N) = L W (N)L 8 (N)L 6 (N)L 4 (N)L 2 (N), 
where the linear differential operators L r (N) are of order r. The first two operators read 

L 2 (N) = Dt 2 + 21 ~ I LH--+ , 1 (3-2) 
2y ' (4-1)4 4i 4(4-1) At 2 K ' 

9 N" 

L 4 (N) = L 4 , - N 2 L^ 2 + — -r, 

lb t 

with Dt = d/dt and: 

4 (24 - 1) o (2414 2 - 2414 + 46) „ 

L 4 = Dt" + 10^ -^-Dt + ^ 5 '-Dt 2 

(t-l)t 2 (4-l) 2 4 2 

(24 - 1) (1224 2 - 1224 + 9) 81 (54 - 1) (5t - 4) 

+ {t- 1) 3 4 3 t+ 16 43 ( t _ !)3 ' 

_ 5 L>4 2 (23 - 324) 9 8 - 174 

4,2 ~ 21 2 " " 2(4-l)4 3 " 8(4- 1)4 4 ' 

The expressions (or forms) of Lq(N), L$(N) and Lio(iV) are given in [25]. The linear differential 
operators F 2n +i(N), which annihilate fjy"^ 1 ^ for n = 0, . . . ,3, are such that: 

F 7 (N) = L 8 (N)L e (N)L 4 (N)L 2 (N), F 5 (N) = L G (N)L A (N)L 2 (N), 

F 3 (N) = L 4 (N)L 2 (N), F^N) = L 2 {N). (3.3) 

Thus we see that the linear differential operator for f^N ^ right divides the linear differential 

operator for /jy^ for n < 3. We conjecture that this property holds for all values of n. We 
thus have a "Russian-doll" (telescopic) structure of these successive linear differential operators. 

3.2 Fuchsian linear differential equations for 

(8) 

The linear differential operator F$(N) (corresponding to f^ N ) has the following factorized form 

F 8 (N) = L 9 (iV)L 7 (iV)L5(7V)L3(iV)Li(iV), 
where the linear differential operators L r (N) are of order r. The first two read: 
Li(N) = Dt, 

r ,™ ~3 A 2t ~ 1 )^ 2 (2 - 154 + 144 2 ) ^ 8t 2_ 15t + 5 / m 1\ 

L 3 (N) = Dt 3 + A) 7v-Dt 2 + ± = '-Dt + a I — k — I — q I N . 

K ' (*-!)< (4-l) 2 4 2 2 (4 -l) 3 4 2 V* 2 t 3 J 

The expressions (or forms) of the linear differential operators L^(N), Lj(N) and Lg(N) are 
given in [25] . 
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Similarly to (|3.3p there is also a Russian-doll (telescopic) structure of these successive linear 
differential operators: 

F 6 (N) = L 7 (N)L 5 (N)L 3 (N)Li(N), F A (N) = L 5 (N)L 3 (N)Li(N), 

F 2 (N) = L 3 (N)Lx{N), F (N) = L^N). (3.4) 

Again, we see that the linear differential operator for /Jv iV ^ right divides the linear differential 
operator for /jy for n < 4. We conjecture that this property holds for all values of n. 

3.3 Direct sum structure 

Not only do the linear differential operators Lj(N) have a factorized Russian-doll structure, 
but we have found that they also have a direct sum decomposition when the integer N is fixed. 
To illustrate this direct sum decomposition, let us write the corresponding linear differential 
operator for fffi N 

F 3 (N) = L A {N)L 2 {N) = M A {N) © L 2 (N), 

where L 2 {N) is the linear differential operator for and where the fourth order opera- 

tor M±(N) is displayed in [25] for successive values of N. One remarks on these successive expres- 
sions that the degree of each polynomial occurring in these linear differential operators M$(N) 
grows linearly with N . 

As a further example consider f( 5 \N,N), where we find that the corresponding linear 
differential operator decomposes as 

F 5 = L 6 (N)L 4 (N)L 2 (N) = M 6 (N) © M 4 (N) © L 2 (N), 

where L 2 (N) is the linear differential operator for fff N , M^(N) is the previous fourth order 
differential operator, and the sixth order operator Mg(iV) has again coefficients whose degrees 
grow with N for successive values of N. There is nothing specific to f§ N and /jy jy : similar 

(n) 

results hold for all the fj^^s, n being even or odd. 

(n) 

In contrast with the Russian-doll way of writing the linear differential operators for fj^ N , 
the direct sum structure, as a consequence of this growing degree, cannot, for generic N, be 
written in a closed form as operators with polynomials coefficients in front of the derivatives. 
This "non-closure" of the direct sum structure will have some consequences when performing 
the scaling limit of these linear differential operators (see Section |4] below) . 

3.4 Equivalence of various Lj(JV)'s and Afj(iV)'s linear differential operators 

We find that the symmetric squar^l of L 2 (N) 

Sym2(L 2 (A0) = Dt* + 3|^W + ^0)-Dt 

1 1-24 N 2 N 2 
-Dt 



2(t-l) 2 t 2 t (t-l)t 2 
and the linear differential operator L 3 (N) are equivalent^ 
L 3 (N)U(N) = V(N)Sym 2 (L 2 (N)) 



5 The symmetric j-th power of a second order linear differential operator having two solutions /i and fi is the 
linear differential operator of order j + 1, which has j '(,..., f(~ /| as solutions. 

6 For the equivalence of linear differential operators, see [1081 11141 1115j . 
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with the following intertwiners: 

U(N) = (t-l) tDt 2 + (St - 1) Dt + 1 + ^— - N 2 , 

V(N) = (t - 1) tDt 2 + (lit - 5) Dt + ^"^"^ - ^T^ 2 - 

Similarly, with the symmetric cube of L 2 (N), we have the equivalence 
L 4 (N)A(N) = B(N)Sjm 3 (L 2 (N)) 

with: 

, ^ , X -k 7 , x n (4U 2 - Alt + 6) 
A(iV) = (t - l)tDt 3 + - (2t - 1) Dt 2 + 1 A[t _ l)t Dt 



9 2t - 1 9(t- 1)N 2 „ 9(2*- 1, , 

+ TVT - T- 1 Dt - - V , 2 7 iV 



2 



(i-l)t 4 £ 8 i 



, , , x o 23 , . 9 21 (6 - 29i + 29i 2 ) 
B(AT) = (t - 1) ^t 3 + — (2t- 1) Dt 2 + -j 1 {t _ 1)t ' Dt 

9 (2t-l)(l25t 2 -125t + 16) _ 9 (i-1) 2 _ 9 (10t-9) 2 
8 (t-l)V 4 t 8 t 2 

More generally, all the L m (iV)'s are (m — l)-symmetric-power of L 2 (N). As a consequen- 
ce their solutions are (m — l)-homogeneous polynomials of the two hypergeometric solutions 
of L 2 (N). 

Similarly, for the linear differential operators occurring in the direct sum, one easily verifies, 
for every integer N, that, for instance, the iW^A^'s are equivalent to the cubic-symmetric-power 
of L 2 (N) 

M 4 (N)Q(N) = S(N)Sym 3 (L 2 (N)) 

where, for N = 0,1, 2: 

Q(0) = (t-l)tDt + t--, 

\3 + 2 



Q(l) = 2 (t - l) 3 t 2 Dt 3 + 3 (3 - 7t + At 2 ) (t - 1) tDt 2 



41 9\ _ 3 2t 2 - 2t + 1 



+ 12t 6 - 28r + — t - - ) Dt + 



2 2; 4 i 

Q(2) = - (i - l) 3 (3 + 8t + 3t 2 ) tDt 3 + - (15 - t - 35t 2 + 15t 3 + 6i 4 ) (i - 1) Dt 2 
3 2 

1 (I8i 5 - 12t 4 - 97i 3 + 577t 2 - 7S8t + 252) 

'-Dt 



24 t 
1 12t 5 + 14t 4 - 260t 3 + 497i 2 - 314t + 24 



16 t 2 

As a further example, one can verify, for every value of the integer N, that the sixth or- 
der operator Mq(N) is equivalent to the fifth symmetric power of L 2 (N). The solutions of 
the linear differential operators M m (N) are also (m — l)-homogeneous polynomials of the two 
hypergeometric solutions of L 2 (N). As a consequence of this direct sum decomposition, the 
solutions (N, N) are (non-homogeneous) polynomials of the two hypergeometric solutions 
of L 2 (N) or, equivalently, f^\ (or the hypergeometric solution of (|3.2p ) and its first derivative. 
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The second order linear differential operator L2(N) is equivalent [27] to the second order linear 
differential operator Le 

L E = AtDt 2 + 4Dt — 

4-1 

corresponding to the complete elliptic integral of the second kind E. As a consequence of the 

(n) 

previously described direct sum decomposition, the fjy N 's can also be written as polynomial 
expressions of the complete elliptic integral of the second kind E and its first derivative E' , or 
alternatively, E and the complete elliptic integral of the first kind K. 

Let us just give here a set of miscellaneous examples of polynomial expressions of various 
form factors. For fffi N , one has 

2/gJ = (K — E) K, 2/g = 1 - 3KE - (t - 2) K 2 , 

6 */2 2 2 = 6* - (2 + 64 2 - lit) K 2 - (154 - 4) KE - 2 (1 + 4) E 2 , 

90i 2 /3,3 = 135* 2 - (1374 3 - 2424 2 + 524 + 8) K 2 

+ (84 3 - 3194 2 + 1124 + 16) KE — 4(1 + 4) (24 2 + 134 + 2) E 2 , 
3 1 504 3 /i 2 4 } = 63004 3 - (324 5 + 25524 2 + 128 + 64404 4 - 111914 3 + 4644) K 2 
+ (1284 5 + 56484 2 - 145194 3 + 10564 + 5764 4 + 256) EK 
- 8 (1 + 4) (164 4 + 584 3 + 3334 2 + 584 + 16) E 2 , 
where E and K are given by (|3.ip . Other examples are given in [25] . 

Miscellaneous remarks. All these remarkable structures are not restricted to diagonal two- 
point correlation functions. We keep on restricting to the isotropic Ising model: for the anisotro- 
pic Ising model one has (for the correlations and may have for the form factors) similar but more 
complicated results involving the complete elliptic integral of the third kind LT (see for instance 
equation (3.35) in H. Au-Yang and J.H.H. Perk [6], or pp. 23-48 in [2], more recently |121j and 
for a sketch of how the algebro-differential structures generalize in that anisotropic case |26j). 

• Further, one can calculate various j-particle contributions fj$ N of the off-diagonal two 
point correlation functions, and verify, again, that they are, in the isotropic case, also polynomial 
expressions of the complete elliptic integrals E and K. For instance: 

C< 2 >(0, l) = 3 -- 1 -(l + s 2 )K- \EK -\{s 2 - 3) (1 + s 2 ) K 2 , 

where s = sinh(2i^). Other miscellaneous examples of such off-diagonal j-particle contributions 

are displayed in [25]. In the anisotropic case polynomial expressions of E and K and complete 

(i) 

elliptic integral of the third kind II could take place for j-particle contributions f^/ N - The 
occurrence of elliptic integral of the third kind and not more involved hyperelliptic integrals is 
still not clear (see after equation (3.20) in [93J the remark on Glasser's unpublished work). This 
work is still in progress. 

• The products of the two-point correlation functions C(N, N) are also solutions of Fuchsian 
linear ODE's. As a consequence the equal-time xx-correlations [98] of the free-fermion zero-field 
XY quantum chain, which are, alternatingly, C(N, N) 2 and C(N, N)C(N+1, N+l), also satisfy 
a Fuchsian linear ODE. 

• Far beyond, recalling Boel, Kasteleyn and Groeneveld papers \20\ [2~T| |4~9] one can see 
that all the two-point correlation functions of Ising models (not necessarily free-fermion Ising 

7 Of course no confusion is possible between the complete elliptic integral K and the usual Ising model tem- 
perature variable sometimes denoted K = J/kT. 
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models!) can be expressed as sums, weighted with ± signs, of products of two-point correlation 
functions. Consequently all the n-point correlation functions of the square Ising model are 
(simple) polynomial expressions of the complete elliptic integrals E and K and, of course, the n- 
point correlation functions of the square Ising model are solutions of Fuchsian linear ODE's. For 
the anisotropic Ising model the n-point correlation functions are solutions of PDE's associated 
with complete elliptic integrals of the third kind (see [26] for a sketch). 

• Recalling the relations (jl.3p , (|1.4p between the two-point correlation functions and the form 
factors we see that, since the isotropic two-point correlation functions and the form factors are 
both polynomial expressions of the complete elliptic integral E and K, relations (jl.3p . (jl.4p 
can be interpreted as an infinite number of quite non trivial identities on the complete elliptic 
integral E and K, for instance: 

oo oo 

C + {N,N) = p[ N > N \ E ,K) = (l-tf'^fNT = (l-t^^Q^^K). 

n=0 n=0 

We have similar identities for the (isotropic) off-diagonal two-point correlations C(M, TV). These 
linear relations on an infinite number of polynomial expressions of the complete elliptic inte- 
grals^) E and K have to be compared with the infinite number of (non-linear) relations on a finite 
number of polynomial identities on the complete elliptic integral E and K which correspond 
to (lA.ip . the quadratic finite difference relations [82|, [84"1 [96| [97] on the two-point correlation 
functions displayed in Appendix lAl 

• At criticality, k = 1, many remarkable and much simpler identities can be obtained, for 
instance the formulg|§ (2.34) in [3] on the next to the diagonal (anisotropic) two-point correlations 
(see also [121] ): 

C(N- 1,N) = C(N,N)cosh(2K)F(l/2,N;N + 1/2, - sinh(2K) 2 ) , 
where F is the hypergeometric function. 



3.5 The elliptic representation of Painleve VI 

The results we have underlined in this section, namely the unexpectedly simple and remarkable 
polynomial expressions for the form factors /Jy V, correspond to the fact that the associated 
linear differential operators are direct sums of operators equivalent to symmetric powers of the 
second order differential operator Le- We already encountered this central key role played 
by the linear differential operator Le, or the hypergeometric second order linear differential 
operator (46) given in [27], in our previous holonomic analysis of the two-point correlation func- 
tions of the Ising model [27] . In order to understand the key role played by Le, or equivalently 
operator L2(N), it is worth recalling (see [78], or for a review [50]) the so-called "elliptic represen- 
tation" of Painleve VI. This elliptic representation of Painleve VI amounts to seeing Painleve VI 
as a "deformation" (see equation (33) in [50]) of the hypergeometric linear differential equation 
associated with the linear differential operator 

C = (1 - t)tDt 2 + (1 - 2t)Dt - -. 

One easily verifies that this linear differential operator is actually equivalent (in the sense of 
the equivalence of differential operators) with Le, or equivalently L2{N). This deep relation 
between elliptic curves and Painleve VI explains the occurrence of Painleve VI on the Ising 

8 In the anisotropic case we could have identities on an infinite number of polynomial expressions of the complete 
elliptic integrals of the first, second and third kind (and hyperelliptic integrals ...?). This work is still in progress. 
9 This formula has first been obtained in [4] at the critical temperature k = 1. 
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model, and on other lattice Yang-Baxter integrable models which are canonically parametrized 
in term of elliptic functions (like the eight-vertex Baxter model, the RSOS models, see for 
instance [15J). One can see in Section 6 of [25], other examples of this deep connection between 
the transcendent solutions of Painleve VI and the theory of elliptic functions, modular curves 
and quasi-modular functions. 

Along this line one should note that other linear differential operators, not straightforwardly 
linked to Le but more generally to the theory of elliptic functions and modular forms (quasi- 
modular forms . . . ) , also emerge in the analysis of the A-extensions of the two-point correlation 
functions of the Ising model, for selected 10 ! values of A: A = cos(7rm/n). This is detailed in 
Appendix [Bj 

4 The scaling limit of ff^ N 

The closed (exact) formulae (|3.3p . (j3.4j) we obtain for the linear differential operators in these 
nested "Russian doll" structures, enable us to take the scaling limit of these linear operators. 
We study this scaling limit in this section and show that the "Russian-doll" structure remains 
valid. The linear differential operators in that "scaled" nested Russian-doll structure remain 
equivalent to the symmetric power of a singled-out second order linear differential operator 
(corresponding to the modified Bessel function). In contrast, in the scaling limit, the direct sum 
of operators decomposition structure is lost, and we explain why. 

in) 

The scaling of the f^^s amounts, on the functions, and on the corresponding differential 
operators, to taking the limit N — > oo and t — > 1, keeping the limit x = N(l — t) finite, or in 
other words, to performing the change of variables t = 1 — x/N, keeping only the leading term 
in N. Performing these straightforward calculations, the linear differential operators in t for 

(n) 

the /^tv's where A?" was a parameter, become linear differential operators in the only scaling 
variable x. 

Calling Fj Cal the scaling limit of the operator Fj(N) we find for j even that 

TTiscal t seal j seal j seal j seal j seal 7-iscal j seal j seal r seal j seal 

1< 8 — ivg 1j 7 Li 5 1j 3 L, x , — 1j 7 iv 5 iv 3 L, x , 

j-iscal t seal j seal t seal 7-iscal j seal j seal jriscal t seal / 4 -1 \ 

^4 = L 5 L 3 L l > b 2 = L 3 L X ) -^0 = L l 1 I 4 - 1 ) 

where (Dx = d/dx) 

L scai = 2x 5 Dx h + l0x 4 Dx* - 2x 3 (7 + 5x 2 ) Dx 3 

+ 2 (-16 + 13x 2 ) x 2 Dx 2 + 2 (5 - 12a; 2 + Ax 4 ) xDx - 10 + 8x 2 - 24x 4 , 
Lf al = 2x 3 Dx 3 + 8x 2 Dx 2 - 2 (x - 1) (x + 1) xDx - 2, 
Lf al = Dx 

and L| cal , L 7 cal are given in [25] . 
Similarly, for j odd, we have 

7-iscal t seal j seal j seal j seal j seal 7-iscal j seal j seal j seal j seal 

-^9 — -^10 ^8 n 6 ^4 Ll 2 1 p 7 — ^8 ^6 ^4 lj 2 ) 

fiscal j seal t seal t seal rascal j seal j seal jpscal rscal /a n\ 

^5 — L S L i L 2 j ^3 — L i L 2 1 *1 — L 2 ) l 4 - 2 J 

where 

L| cal = 16x 4 £>x 4 + 9Qx 3 Dx 3 + 40 (2 - x 2 ) x 2 Dx 2 + 8 (x 2 - 2) xDx + 9x 4 - 8x 2 + 16, 
L| cal = AxDx 2 + 4Dx - x 

and L^ 1 , L| cal , Lg cal are given in [25] . 



3 For generic values of A, the A-extension C(M, N; A) are not holonomic. 
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Thus, we see that the scaled operators i^ cal have a "Russian-doll" structure straightforwardly 
inherited from the one for the lattice operators Fj(N). 

Consider the linear second-order differential operator corresponding to the modified Bessel 
function K n (x/2) for n = 0, namely: 

„ „ 9 Dx 1 

B = Dx 2 + -. (4.3) 

x 4 

We recognize, in this linear differential operator, the exact identification with the scaled diffe- 
rential operator Ff cal = L| cal . We find that the symmetric square of the linear differential 
operator B, and the scaled operator L| cal are equivalent: 

L s 3 cal (xDx 2 + 2Dx - x) = (2x i Dx 2 + l2x 3 Dx - 2x A + 8x 2 )Sym 2 (5). 

Similarly, the symmetric third power of the linear differential operator B, and the scaled oper- 
ator L| cal are equivalent, and, more generally, the symmetric j-th. power of (|4.3|) and the scaled 
operator LJ^\ are equivalent: 

Lf*\ ~ Sjm j (B). 

Recall that the linear differential operators Fj(N), corresponding to the form factors fjp N , 

can be written as direct sums only when the integer iV is fixed. At the scaling limit, this feature 

disappears for the scaled linear differential operators i^? cal which have no direct sums. Therefore 

while the scaling limit preserves the Russian-doll (telescopic) structure (see (|3.3[) . (|4.2p ) and 

also preserves the fact that the various operators in this Russian-doll (telescopic) structure are 

equivalent to symmetric powers of an operator (|4.3p which replaces the operator Le, the direct 

(i) 

sum structure is lost. As a consequence the scaling of the f^j^'s cannot be seen as simple 
polynomials of modified Bessel functions. 

There is one exception that concerns fjy\i- Its scaled linear differential operator _F| cal , has 
the non shared property of being equivalent to the direct sum of Dx with the symmetric square 
of (|4.3j) . namely: 

F scal = L scal Q L scal ^ Dx Q Sym 2 (£). 

(2) 

From this equivalence, one immediately deduces the expression of the scaling of the f^f N as 
a quadratic expression of the modified Bessel functions of x/2 which actually identifies with 
formula (2.31b)-(3.151) in [123] . 

The occurrence of modified Bessel functions, emerging from a confluence of two regular 
singularities of the complete elliptic integrals E and K, or from the hyper geometric function 
2-F1, should not be considered as a surprise if one recalls the following limit of the hyper geometric 
function 2 Fx yielding confluent hypergeometric functions \F\. These confluent hypergeometric 
functions, %Fi, are nothing but modified Bessel functions [37] 

2F1 fa,p,b;^\ — > iF 1 (a,b;z) when p —> 00, 

z u f 1 

U v z) = F x \v + -, 2v + 1; 2z 

Remark. It was shown, in Section [31 as a consequence of the decomposition of their linear 
differential operators in direct sums of operators equivalent to symmetric powers of Le, that 
the functions /jyL are polynomial expressions of E and K functions. Therefore their singularities 
are only the three regular points t = 0, t = 1 and t = 00. The scaling limit (t = 1 — x/N, t — > 1, 



14 



S. Boukraa, S. Hassani, J.-M. Maillard and N. Zenine 



N — > oo) corresponds to the confluence of the two regular singularities t = and t = oo, yielding 
the, now, irregular singularity at x = oo. The occurrence of irregular singularities with their 
Stokes phenomenon, and, especially, the loss of a remarkable direct sum structure, shows that 
the scaling limit is a quite non-trivial limit. 

Contrary to the common wisdom, the scaling limit does not correspond to more "fundamen- 
tal" symmetries and structures (more universal . . . ) : this limit actually destroys most of the 
remarkable structures and symmetries of the lattice model^l. 



5 Bridging with other formula of form factors 
in the scaling limit: work in progress 

The Ising Form Factors in the scaling limit as they can be found in Wu, McCoy, Tracy, 
Barouch [123j read 

-i poo poo 

and: 

i poo poo 2k+l —ty-j ^fc k 

^ - i-v^L *™ n j^iij^ms, - d- 

The Ising Form Factors in the scaling limit are also given, in many field theory papers, as follows 
(yj = coshOj, see (9) in (7j, see also Mussardo [H6]): 

^ = ¥ /_! • • • /.I n n ^ (^) ■ (-> 

It remains to show that those expressions are actually solutions of the scaled linear differential 
operators displayed in the previous Section HI namely (|4.ip and (|4.2p . A straight check yields 
too large formal calculations. Our strategy should rather be to obtain the series expansions of 
the n-fold integrals (|5.ip or (15. 2j) . in the t variable for (|5.ip . or the r variable for (|5.2p . and 
check that these series expansions are actually solutions of the non-Fuchsian linear differential 
operators (|4~T1) and (|4~2l) of Section [U 

These checks will show that the expressions f)5. 1[) . ()5.2[) of the n-fold integrals of the scaled 
form factors are actually solutions of linear differential operators with an irregular singularity at 
infinity and with a remarkable Russian-doll structure but no direct sum structure. This will indi- 
cate that such expressions (|5.ip . (|5.2p generalise modified Bessel functions but cannot be simply 
expressed in terms of polynomial expressions of modified Bessel functions. The interpretation 
of these expressions in terms of r- functions (Hirota equations, hierarchies, . . . ) and the link 
between these Russian-doll structures and Backhand transformations or Hirota transformations 
remains to be done in detail. 



6 Other n-fold integrals: from diagonal correlation functions 
to the susceptibility of the Ising model 

The study of two-point correlation functions (even n-points correlations . . . ) can be seen as 
a "warm-up" for the truly challenging problem of the study of the susceptibility of the Ising 

n This kind of results should not be a surprise for the people working on integrable lattice models, discrete 
dynamical systems, or on Painleve equations [85lll06| . 
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model and its associated n-fold integrals, the x^ ( see next Section [7] below) . Staying close to 
the diagonal correlation functions we have introduced a simplification of the susceptibility of the 
Ising model by considering a magnetic field restricted to one diagonal of the square lattice [28] . 
For this "diagonal susceptibility" model [28], we benefited from the form factor decomposition 
of the diagonal two-point correlations C(N,N), that has been recently presented [25], and 
subsequently proved by Lyberg and McCoy [74]. The corresponding n-fold integrals x* were 
found to exhibit remarkable direct sum structures inherited from the direct sum structures of 
the form factor [25^ I28j. The linear differential operators of the form factor [25] being closely 
linked to the second order differential operator Le (resp. Lk) of the complete elliptic integrals E 
(resp. K), this "diagonal susceptibility" |28j is also closely linked to the elliptic curves of the 
two-dimensional Ising model. By way of contrast, we note that the singularities of the linear 
ODE's for these n-fold integrals [28] x^ are quite elementary (consisting of only n-th roots 
of unity) in comparison with the singularities we will encounter below with the quite simple 
integrals (|7T4j) . 

Using the form factor expansions (jl.5p and (|1.6p , the A-extension of this diagonal susceptibility 
can be written as 

oo 

k B T Xd ±{t-A) = {i-t) 1 l i Y. x ^% 

j 

where the sum is over j even (resp. odd) for T below (resp. above) T c and where 

oo 

X<2± = /aT,Jv(*)- 
N= 



By use of the explicit expressions (|2.5|) . and f|2.6[) . for fjp N we find explicitly, for T < T c , that 



_(2n) = _^_J_ f 1 / ,1 TT J l+* n Sl---S2n -TT ( X 2j -l(l ~ X 2j )(l - tx 2j ) \ 1/2 

" (n!)2 J "'J 11 Xk l-t"x 1 ---x 2n f} 1 \x 2j (l-x 2j - 1 )(l-tx 2j - 1 )J 



X _ (1 - tx 2 j-\X 2k ) 2 (X2j~l - X2k-l) 2 (x2j - X 2 k) 2 

l<j<nl<k<n l<j<fc<n 



and, for T > T c , that: 



t n( n +i)) ^ r 1 r 1 2n+1 
^ 2;+1 = n!(n TTy.^Jo "'J 

i I j.n+1/2^. _ n + 1 

L + t ' Xl ■ ■ ■ X 2n+ l yj x/ 2 



fe=l 

n+1 
n+1 

x JJ((l-o;2i-i)(l-te2i-i)ac2i-ir 1/2 II II ^ " te 2i-1^2fe)" 2 
j=l l<i<n+l l<fc<n 

X {X2j-1 ~ X 2k -lf {X2j-X 2 kf- 

l<j'<fc<n+l l<j<k<n 

We have also found [28], for j = 1, ... ,4, that the Xd± s satisfy Fuchsian linear differential 
equations which have a Russian-doll nesting just as was found for the Xn's in [126^11271 11281 [T29] . 
In the case of these j-particle components of the "diagonal" susceptibility, we can see that 
this Russian-doll nesting of the corresponding linear differential operators is straightforwardly 
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inherited [25], not from the Russian-doll nesting of the diagonal form factors /^(iV, iV)'s (this 
is not sufficient), but from their direct sum (of operators equivalent to symmetric powers) 
decomposition. 

Direct sum decompositions of the xjj_ ana - Xd+ +1 ^ are straightforwardly inherited [2H] from 
direct sums of the /Jrf N thus yielding a scenario for the direct sum decompositions of the 
"true" ■ However, recalling the non-unicity of the {G^n , F^fy (see Section [2]), the 
direct sum decomposition of the x^ can De seen as a canonical one, when the direct sum 
decomposition of the ^( 2n+1 ) is not. 



7 Other n-fold integrals linked to the susceptibility 
of the Ising model 

The susceptibility x °f the square lattice Ising model has been shown by Wu, McCoy, Tracy and 
Barouch [123] to be expressible as an infinite sum of holomorphic functions, given as multiple 
integrals, denoted x (n) , that is kTx = E X (n) - B. Nickel found (87] EE] that each of these x (n) 's 
is actually singular on a set of points located on the unit circle |s| = |sinh(2X)| = 1, where 
K = J/kT is the usual Ising model temperature variable. 

These singularities are located at solution points of the following equations 

L = 2 ( s+ 1 -) =u k + ±+u™ + ± (7.1) 
with u 2n+l = 1, — n < m, k < n. 



From now on, we will call these singularities of the "Nickelian type", or simply "Nickelian 
singularities" . The accumulation of this infinite set of singularities of the higher-particle com- 
ponents of x( s ) on the unit circle \s\ = 1, leads, in the absence of mutual cancelation, to 
some consequences regarding the non holonomic (non D-finite) character of the susceptibility, 
possibly building a natural boundary for the total x( s )- However, it should be noted that 
new singularities, that are not of the "Nickelian type", were discovered as singularities of 
the Fuchsian linear differential equation associated [1261 1128^ 1129] with x^ an d as singu- 
larities of x^ itself [29] but seen as a function of s. They correspond to the quadratic 
polynomial 1 + 3w + Aw 2 where 2w = s/(l + s 2 ). In contrast with this situation, the Fuch- 
sian linear differential equation, associated [127| with x^\ does not provide any new singula- 
rities. 

Some remarkable "Russian-doll" structure, as well as direct sum decompositions, were found 
for the corresponding linear differential operators for x^ an d X ■ I n order to understand 
the "true nature" of the susceptibility of the square lattice Ising model, it is of fundamen- 
tal importance to have a better understanding of the singularity structure of the n-particle 
contributions x > an d also of the mathematical structures associated with these x^ n \ namely 
the infinite set of (probably Fuchsian) linear differential equations associated with this infinite 
set of holonomic functions. Finding more Fuchsian linear differential equations having the 
%( n )'s as solutions, beyond those already found [1261 1127] for x^ an d X^\ probably requires 
the performance of a large set of analytical, mathematical and computer programming "tours- 
de-force" . 

As an alternative, and in order to bypass this "temporary" obstruction, we have developed, 
in parallel, a new strategy. We have introduced [29] some single (or multiple) "model" integrals 
as an "ersatz" for the x^'s as far as the locus of the singularities is concerned. The x^' s are 
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defined by (n — l)-dimensional integrals |88[ [94"] I124j (omitting the prefactoJ^l) 



where 



n\ - LJ - In 2tt 

j=i J0 \j=i 



(n-l)/2 



(n ) n 2s T;^ )/2) . «w=^ 

\i=l / l<«<j<n 8 j 1 - [I ^ 



i=l 



with 



2w 



(7.3) 



1 - 2w cos(0i) + y/ (1 - 2w cos((pi)y - Aw 2 
1 



^(l-2wcos(^)) 2 -4u; 2 ' j=i 



The two families of integrals we have considered in [29] are very rough approximations of the 
integrals (|7.2h . For the first family^!, we considered the n-fold integrals corresponding to the 
product of (the squar^l of the) y^'s, integrated over the whole domain of integration of the 4>i 
(thus getting rid of the factors and R^). Here, we found a subset of singularities occurring 
in the x as we ^ as the quadratic polynomial condition 1 + 3w + Aw 2 = 0. 

For the second family, we discarded the factor and the product of y^'s, and we restricted 
the domain of integration to the principal diagonal of the angles <f>i {<pi = (fa = • ■ • = 4> n -i)- 
These simple integrals (over a single variable), were denoted [23] 



1 2 



2k 



D n\ nlj 2^rl-x n - 1 (<p)x({n-l)<P) , K ' J 

where x(4>) is given by (|7.3p . 

Remarkably these very simple integrals both reproduce all the singularities, discussed by 
Nickel [871 [88], as well as the quadratic roots of 1 + 3w + Aw 2 = found [HH [J29] for the 
linear ODE of ■ O ne should however note that, in contrast with the x^ n \ no Russian-doll, or 
direct sum decomposition structure, is found for the linear differential operators corresponding 

(n) 

to these simpler integrals & D ■ 

We return to the integrals (|7.2p where, this time, the natural next step is to consider the 
following family of n-fold integrals 



71 

n— 1 „2-7r ii ( n \ 1 ll "^i 

i=i 



12 The prefactor reads (1 — s 4 ) 1/4 /s for T > T c and (1 — s 4 ) 1 / 4 for T <T C and in terms of the w variable. 
13 Denoted Y {n) {w) in [29]. 

14 Surprisingly the integrand with (11^=1 Vj) 2 yields second order linear differential equations [29], and conse- 
quently, we have been able to totally decipher the corresponding singularity structure. By way of contrast the 
integrand with the simple product (I~[j=i Hi) yields linear differential equations of higher order, but with identical 
singularities [29] . 
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which amounts to getting rid of the (fermionic) factor (G^) 2 in the n-fold integral (|7.2p . This 
family is as close as possible to (|7.2|) . for which we know that finding the corresponding linear 
differential ODE's is a huge task. The idea here is that the methods and techniques we have 
developed [1261 1129j for series expansions calculations of and x^ > seem to indicate that the 
quite involved fermionic term (G^) 2 in the integrand of (|7.2p should not impact "too much" 
on the location of singularities of these n-fold integrals (17. 2j) . This is the best simplification 
of the integrand of (|7.2p for which we can expect to retain much exact information about the 
location of the singularities of the original Ising problem. However, we certainly do not expect 
to recover from the n-fold integrals (|7.5p the local singular behavior (exponents, amplitudes of 
singularities, etc . . . ). Getting rid of the (fermionic) factor (G^) 2 are we moving away from the 
elliptic curves of the two-dimensional Ising model? Could it be possible that we lose the strong 
(Russian-doll, direct sum decomposition) algebro-differential structures of the corresponding 
linear differential operators inherited from the second order differential operator Le (resp. Lk) 
of the complete elliptic integrals E (resp. K), but keep some characterization of elliptic curves 
through more "primitive" (universal) features of these n-fold integral like the location of their 
singularities? 

In the sequel, we give the expressions of and the Fuchsian linear differential equa- 

tions for for n = 3 and n = 4. For n = 5,6, the computation (linear ODE search of 
a series) becomes much harder. Consequently we use a modulo prime method to obtain the 
form of the corresponding linear ODE with totally explicit singularity structure. These results 
provide a large set of "candidate singularities" for the x ■ From the resolution of the Landau 
conditions |29[ [35] for (|7.5p , we have shown that the singularities of (the linear ODEs of) these 
multiple integrals actually reduce to the concatenation of the singularities of (the linear ODEs 
of) a set of one-dimensional integrals. We discuss the mathematical, as well as physical, inter- 
pretation of these new singularities. In particular we can see that they correspond to pinched 
Landau-like singularities as previously noticed by Nickel |89j . Among all these polynomial sin- 
gularities, the quadratic numbers 1 + 3w + Aw 2 = are highly selected. We will show that 
these selected quadratic numbers are related to complex multiplication for the elliptic curves 
parameterizing the square Ising model. 

(n) 

We present the multidimensional integrals & H and the singularities of the corresponding 
linear ODE for n = 3, . . . , 6, that we compare with the singularities obtained from the Landau 
conditions. We have shown [30] that the set of singularities associated with the ODEs of the 
multiple integrals & H reduce to the singularities of the ODEs associated with & finite number of 
one- dimensional integrals. Section [9] deals with the complex multiplication for the elliptic curves 
related to the singularities given by the zeros of the quadratic polynomial 1 + 3w + Aw 2 . 

8 The singularities of the linear ODE for <1>^ 



For the first two values of n, one obtains 



(n) 

For n > 3, the series coefficients of the multiple integrals $W are obtained by expanding in 
the variables Xi and performing the integration (see Appendix A of [30]). One obtains 

oo oo 

= n\ £ " M(2 " 6 pfi )w n ^a n (k,p), (8.1) 

' fe=0p=0 



where a(k,p) is a 4F3 hypergeometric series dependent on w. 
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The advantage of using these simplified integrals (|7.5p instead of the original ones (|7.2p is 
twofold. Using (|8.1|) the series generation is straightforward compared to the complexity related 
to the • As an illustration note that on a desk computer, $^ are generated up to u; 200 in less 
than 10 seconds CPU time for all values of n, while the simplest case of the x^ n \ namely x^ 3 \ 
took three minutes to generate the series up to w 200 . This difference between the <3?^ and x^ 
increases rapidly with increasing n and increasing number of generated terms. We note that 
for the <3?^ quantities and for a fixed order, the CPU time is decreasing^ with increasing n. 
For x^- the opposite is the case. The second point is that, for a given n, the linear ODE can 
be found with less terms in the series compared to the linear ODE for the x^ ■ Indeed for x^ 3 \ 
360 terms were needed while 150 terms were enough for ■ The same feature holds for x^ 
and <J>^ (185 terms for x^ an d 56 termj^l for 

With the fully integrated sum (|8.ip . a sufficient number of terms is generated to obtain 
the linear differential equations. We succeeded in obtaining the linear differential equations, 
respectively of minimal order five and six, corresponding to 3>S^ and These linear ODE's 
are given in Appendix ICl 

For $y (n > 5), the calculations, in order to get the linear ODEs become really hugj^l. For 
this reason, we introduce a modular strategy which amounts to generating long series modulo a 
prime and then deducing the ODE modulo that prime. Note that the ODE of minimal order is 
not necessarily the simplest one as far as the required number of terms in the series expansion 
to find the linear ODE is concerned. We have already encountered such a situation \127\ 128]. 
For <J>^ (resp. < &^), the linear ODE of minimal order is of order 17 (resp. 27) and needs 8471 
(resp. 9272) terms in the series expansion to be found. 

Actually, for (resp. we have found the corresponding linear ODEs of order 28 

(resp. 42) with only 2208 (resp. 1838) terms from which we have deduced the minimal ones. The 
form of these two minimal order linear ODEs obtained modulo primes is sketched in AppendixO 
In particular, the singularities (given by the roots of the head polynomial in front of the highest 
order derivative), are given with the corresponding multiplicity in AppendixO Some details 
about the linear ODE search are also given in Appendix [Cl 

(7) 

We have also obtained very long series (40000 coefficients) modulo primes for but, un- 
fortunately, this has not been sufficient to identify the linear ODE (mod. prime) up to order 100. 

In) 

The singularities of the linear ODE for the first $ H are respectively zeros of the following 
polynomials (besides w = oo): 

n = 3, w (1 - 16w 2 ) (1 - w) (1 + 2w) (l + 3w + Aw 2 ) , 
n = 4, w (1 - 16w 2 ) (1 - 4u> 2 ) , 

n = 5, w (1 - 16w 2 ) (l - w 2 ) (1 + 2w) (l + 3w + Aw 2 ) (l - 3w + w 2 ) 
x (l + 2w- Aw 2 ) (1 + Aw + 8w 2 ) (l-7w + 5w 2 - Aw 3 ) 

x (1 - w - 3w 2 + Aw 3 ) (1 + 8w + 20w 2 + 15w 3 + Aw A ) , (8.2) 
n = 6, w (1 - 16u; 2 ) (l - Aw 2 ) (l - w 2 ) (l - 25-u; 2 ) (l + 3w + Aw 2 ) 

x (1 - 9-u; 2 ) (l - 3w + 4-u; 2 ) (l - Ww 2 + 29w A ) . (8.3) 

For n = 7 and n = 8, besides modulo primes series calculations mentioned above, we also 
generated very long series from which we obtained in floating point form, the polynomials given 
in Appendix [D] (using generalised differential Pade methods). 

15 This can be seen from the series expansion (|8.1[) . Denoting Ro the fixed order, one has n(p + k) < Ro, while 
the CPU time for the series generation of a n (k,p) is not strongly dependent on n. 

16 From now on, for even n, the number of terms stands for the number of terms in the variable x = w 2 . 
17 Except the generation of long series which remains reasonable. 
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(n) 

If we compare the singularities for the ODEs for the <& H to those obtained with the "Diagonal 
model"(lf| presented in [29], i.e. for the ODEs for the one sees that the singularities of the 

(n) 

linear ODE for the "Diagonal model" are identical to those of the linear ODE of the $ H for 

n = 3, 4 (and are a proper subset to those of <J>^ for n = 5, 6). The additional singularities for 
n = 5, 6 are zeros of the polynomials: 

n = 5, (l + 3w + Aw 2 ) (l + Aw + 8w 2 ) (l - 7w + 5w 2 - Aw 3 ) , 
n = 6, (l + 3w + Aw 2 ) (1 - 3w + Aw 2 ) (l - 25w 2 ) . 

For ra = 7, the zeros of the following polynomials (among others) are singularities which are not 
of Nickel's type (|7.ip and do not occur for $^ } : 

1 + 8w + 15w 2 - 21w' 3 - 60u> 4 + 16w 5 + 96w 6 + 6Aw 7 , 
1-Aw- 16w 2 - A8w 3 + 32u> 4 - 128w 5 . 

The linear ODEs of the multiple integrals <fr^ thus display additional singularities for n = 5, 6 
and n = 7 (n = 8 see below) compared to the linear ODE of the single integrals . 

(n) 

We found it remarkable that the linear ODEs for the integrals & D display all the "Nickelian 

singularities" (17. ip , as well as the new quadratic numbers 1 + 3w + Aw 2 = found for ■ It is 

M (n) 
thus interesting to see how the singularities for & D are included in the singularities for & H and 

whether the new (with respect to singularities can be given by one-dimensional integrals 

similar to 3>r, . Let us mention that the singularities of the linear ODE for $^ (resp. 3>S ) are 

remarkably also singularities of the linear ODE for <E>^ (resp. $^)- In [30] it was shown how 
this comes about and how it generalizes. For this, we had to solve the Landau conditions [30] 
for the n-fold integrals (|7.5p . 

9 Bridging physics and mathematics 

In a set of papers [26\ 127] and in the previous sections, we have underlined the central role played 
by the elliptic parametrization of the Ising model, in particular the role played by the second 
order linear differential operator Le (or Lk) corresponding to the complete elliptic integral E 
(or K), and the occurrence of an infinite number of modular curves [2S], canonically associated 
with elliptic curves. We are getting close to identify the lattice Ising model, (or more generally 
Baxter model), with the theory of elliptic curves. In such an identification framework one may 
seek for "special values" of the modulus k that could have a "physical meaning", as well as 
a "mathematical interpretation" (beyond just being singularities), as singularities of the x ■ 

9.1 Revisiting the theory of elliptic curves with a physics viewpoint 

The deep link between the theory of elliptic curves and the theory of modular forms is now 
well established [105J . More simply the crucial role of the modular group in analysing elliptic 
curves is well known. For that reason seeking "special values |71j " of the modulus k, that might 
have a "physical meaning" as well as a mathematical meaning, as singularities of the X \ h 
may be interesting to, alternatively, introduce the modular function called the j-function which 
corresponds to Klein's absolute invar iant multiplied by (12) 3 = 1728 



3 Not to be confused with the "diagonal susceptibility" and the corresponding [28] n-fold integrals xT 
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and, alternatively, seek for "special values" of the j-function (|9.ip . since it automatically takes 
into account the modular symmetry group of the problem. The modular group requires one to 
introduce the period ratio and the nome of the elliptic functions. The elliptic nome, defined in 
terms of the periods of the elliptic functions, reads 

/ K{l-k 2 )\ . . 

q = exp I — 7r K ^ 2 j J = exp(ivrr), (9.2) 



where r is the half period ratio 19 !. 

The SL(2, Z) transformations of the modular group 

ar + b 



ct + d 



which preserve the j-function (|9.ip . should not be confused with isogenics of elliptic curves like 
the Gauss or Landen transformations 

r - 2r, r - I (9.3) 
and, more generally (n integer) 

7" 

r — ► nr, r — > — . (9.4) 

n 

which actually modify the j-function (|9.ip . but are "compatible" with the lattice of periods (the 
inclusion of one lattice into the other one) . 

Roughly speaking, and as far as the elliptic curves of the Ising model (resp. Baxter model) 
are concerned, the SL(2, Z) transformations of the modular group are invariance symmetries 
(reparametrizations), while the transformations (|9.4h are highly non-trivial covariants that we 
will see as exact representations of the renormalization group. 

9.2 Landen and Gauss transformations as generators 
of the exact renormalization group 

Let us consider the complete elliptic integral K(k) defined as: 

K(k) = 2-F1 (1/2, 1/2; 1; k) . 

Two relations between K(k), evaluated at two different modulus, can be found in, e.g. [36] and 
rear" 20 '- 



(l+J0/f(t 2 )=A'( Ir ^ ? ). (9.6) 

The arguments in K, on the right-hand-side of (|9.5p . (|9.6p . are the square of the modulus k 
transformed by the so-called (descending) Landen or (ascending) Landen (or Gauss) transfor- 
mations: 

h , h , WjEg (97 ) 



19 In the theory of modular forms q 2 is also used instead of q. In number theoretical literature the half-period 
ratio is taken as —it. 

20 Note that the relation (|0 . 5 f) is valid for any value of the modulus k, while the validity of (|9.6|l is restricted to 



|fc| < 1. 
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1 = T+k' (98) 

A sequence of such transformations can be used to evaluate (numerically) , in a rapidly convergent 
way, the elliptic integrals from iterations of (19, 7p or of (|9.8j) . Changing k to the complementary 
modulus k! = \f\ — k 2 , and likewise for the transformed k, the half period ratio transforms 
through (|22D, dSSD, like (l9~3j) . 

The real fixed points of the transformations (19. 7h and (19. 8p are A; = (the trivial infinite 
or zero temperature points) and fc = 1 (the ferromagnetic and antiferromagnetic critical point 
of the square Ising model). Iterating (|9.7p or (|9.8p . one converges, respectively to A; = or 
A; = 1. In terms of the half period ratio, this reads, respectively, r = oo and r = which 
correspond to a degeneration of the elliptic parametrization into a rational parametrization. In 
view of these fixed points, it is natural to identify the transformations (|9.7p or (|9.8p . and, more 
generally, any transformational r — > iit or t — > r/n (n integer), as exact generators of the 
renormalization group. It is a straightforward exercise, using the identities (|9.5p . (19. 6p . to write 
a "renormalization recursion" on the internal energy U of the Ising model 

where c and s denote cosh(2i^T) and sinh(2i^) respectively. 



9.3 Complex multiplication of elliptic curves 

and fixed points of Landen transformations 

Since we are interested in singularities in the complex plane of some "well-suited" variable (s, k, 
w), one should not restrict (|9.7p and (|9.8p as transformations on real variables, restricting to 
real fixed points of these transformations, but actually consider the fixed points of these trans- 
formations seen as transformations on complex variables. 

For instance, if one considers (19. 8p as an algebraic transformation of the complex variable k 
and solve k\ — k 2 = 0, one obtains: 

£;(l-£;)(£; 2 + 3A; + 4) = 0. 

The roots of 

k 2 + 3k + i = 0, (9.9) 

are (up to a sign) fixed points of (|9.8p . We thus see the occurrence of additional non-trivial com- 
plex selected values of the modulus k, beyond the well-known values k = 1,0, oo, corresponding 
to degeneration of the elliptic curve into a rational curve, and physically, to the critical Ising 
model and to (high-low temperature) trivializations of the model. 

Of course, when extending (|9,8p to complex values, one can be concerned about keeping track 
of the sign of k\ in (|9.8p in front of the square root \fk. Reference |30] provides a similar fixed 
point calculation for (|9.8p extended to complex values, but for a representation of (|9.8p in term 
of the modular j'-function. Such calculations single out the remarkable integer value j = (— 15) 3 , 
which is known to be one of the nine Heegner numbers (see [30]) . It is important to note that 
this representation of (|9.8p in term of the modular j-function is the well-known fundamental 
modular curve symmetric in j and ji (see [301 [5T| [75] ) 

J 2 Ji 2 - (j + ii)(j' 2 + 1487JJ! + ft) - 12 • 30 6 (j + h) 



See relation (1.3) in [111] . 
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+ 3 • 15%16f - 4027jji + 16jf) + 8 • 30 9 = 0. 

which represents, at the same time, the Landen and Gauss transformations (|9.3p as a consequence 
of the modular invariance (r <-> 1/r). 

A straightforward calculation of the elliptic nome (|9.2p gives, for the polynomial (|9.9|) and 
the polynomial deduced from the Kramers- Wannier duality k — > 1/k respectively, exact values 
for r, the half period ratio, as very simple quadratic numbers: 

±3 + iV7 ±l + iy/f 
ri = 4 , r 2 = . 

These quadratic numbers correspond to complex multiplication and to j = (— 15) 3 . These two 
quadratic numbers are such that 2t\ =f 1 = r 2 . Let us focus on r 2 for which we can we write 

2 

r = 1 - -. 
r 

Taking into account the two modular group involutions r — ► 1 — r and r — > 1/r, we find that 
1 — 2/r is, up to the modular group, equivalent to r/2. The quadratic relation r 2 — r + 2 = thus 
amounts to looking at the fixed points of the Landen transformation r — ► 2r up to the modular 
group. This is, in fact a quite general statement: the complex multiplication values can all be 
seen as fixed points, up to the modular group, of the generalizations of Landen transformation, 
namely r — > nr for n integer (here ~ denotes the equivalence up to the modular group): 

2 n 1 + iy/4n - 1 

t — t + n = U or r=l ~ nr, r = . 

r 2 

Complex multiplication corresponds to integer values of the modular j'-function (as in the case 
of the Heegner numbers sec |30j). 

For elliptic curves in field of characteristic zero, the only well-known selected set of values 
for k corresponds (besides k = 0,1, oo) to the values for which the elliptic curve has complex 
multiplication [77], and we see these selected values, here, as fixed points, in the complex plane, 
of transformations (isogenies) that are exact representations of generators of the renormalization 
group. 

It is now totally natural to see if the singularities we have obtained for the n-fold inte- 
grals (|7.5j) , can be interpreted in the framework of elliptic curve theory, in terms of this physically, 
and mathematically, highly selected set of values for elliptic curves, namely complex multiplica- 
tion values. 



9.4 Complex multiplication for 1 + 3w + 4w 2 = 

Let us consider the first unexpected sing ularities 1 + 3w + 4w 2 = we found [T26| [T29] for the 
Fuchsian linear differential equation of x \ an d also found in other n-fold integrals of the Ising 
class [29J. This polynomial condition reads in the s variable, (2s 2 + s + l) (s 2 + s + 2J =0. We 
have shown [29] that itself is not singular at the roots of the first polynomial whose roots 
are such that \s\ < 1, but is actually singular at the roots of the second polynomial. In the 
variable k = s 2 , these singularities read: 

(4k 2 + 3k + l){k 2 + 3/c + 4) = 0. (9.10) 

The second polynomial has actually been seen to correspond to fixed points of the Landen 
transformation (see (|9.9p ). Note that the two polynomials in (|9.10p are related by the Kramers- 
Wannier duality k — > 1/k (and therefore both correspond to the same value of the modular 
j-function: j = (— 15) 3 ). 
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In other words we see that the selected values 1 + 3w + Aw 2 = 0, occurring in the (high- 
temperature) susceptibility of the Ising model as singularities of the three-particle term x^ 3 \ 
actually correspond to the occurrence of complex multiplication on the elliptic curves of the 
Ising model, and can also be seen as fixed points of the renormalization group when extended to 
complex values of the modulus k. 

Let us note that the occurrence of Heegner numbers and complex multiplication has already 
occurred in other contexts, even if the statement was not explicit. In the framework of the 
construction of Liouville field theory, Gervais and Neveu suggested [46] new classes of critical 
statistical models (see Appendix[E|), where, besides the well-known iV-th root of unity situation, 
they found the following selected values of the multiplicative crossing |103j t: 

t = e WlW3)/2 = ie -nV3/2 } (9 n) 

t = e i7T{1+i) = -e - \ (9.12) 

If one wants to see this multiplicative crossing [30l [52l E3 EO] as a modular nome (see [30]), the 
two previous situations actually correspond to selected values of the modular j-function namely 
j((l + iV3)/2) = (0) 3 for (I9TTTD . and j(l + i) = (12) 3 for (I9J2D . which actually correspond 
to Heegner numbers and complex multiplication [77]. It is however important not to feed the 
confusion already too prevalent in the literature, between a "temperature-like" nome like (|9.2p 
and a multiplicative crossing modular nome (see Appendix [E]). In the Baxter model [IT], the 
first is denoted by q and the second one by x. In fact one probably has, not one, but two 
modular groups taking place, one acting on the "temperature-like" nome q and the other one 
acting on the multiplicative crossing x. We will not go further along this quite speculative 
line which amounts to introducing elliptic quantu m g roups [79] and (see Appendix [E]) elliptic 
gamma functions (generalization of theta functions 22 !) which can be seen (39] as "automorphic 
forms of degree 1" , when the Jacobi modular forms are "automorphic forms of degree 0" and 
are associated (up to simple semi-direct products) to SL(3, Z) instead of SL(2, Z). 



9.5 Beyond 1 + 3w + 4w 2 = 

As a consequence of the fact that the modular j-function is a function of w 2 , the quadratic 
polynomial condition 1 — 3w + Aw 2 = 0, corresponds to the same selected values of the modular 
j-function as 1 + 3w + Aw 2 = 0, namely j = (— 15) 3 . The quadratic polynomial 1 — 3w + Aw 2 = 
actually occurs in the singularities of the linear ODE for <3?^ (and all the higher < 3?^ 2n \ if one 
believes formulas (28) and (29) in [30] ) . 

In view of the remarkable mathematical (and physical) interpretation of the quadratic values 
l+3w+Aw 2 = 0, (and also 1 — 3w+4w 2 = 0) in terms of complex multiplication, or fixed points of 
the renormalization group, it is natural to see if such a "complex multiplication" interpretation 
also exists for other singularities of x^ n \ an d as a first step, for the singularities of the linear 
differential equations of our n-fold integrals (|7.5p . that we expect to identify, or at least, have 
some overlap with the singularities of the x • 

We have found two other polynomial conditions which correspond to remarkable integer values 
of the modular j-function. The singularities 1 — 8w 2 = correspond to j = (12) 3 and r = ±1 +i 
(see |30j). They correspond to "Nickelian singularities" ()7.ip for x^ (and thus $S ) and to "non- 
Nickelian singularities" for and $^ ^ Another polynomial condition is 1 — 32m; 2 = 0, which 
give "non-Nickelian singularities" that begin to appear at n = 10. These singularities correspond 
to the integer value of the modular j-function, j = (66) 3 and to r = 2i and r = —4/5 + i2/5. 

22 The partition function of the Baxter model can be seen as a ratio and product of elliptic gamma functions 
and theta functions, it is thus naturally expressed as a double infinite product. Similar double, and even triple, 
products appear in correlation functions of the eight vertex model [59j [61] . 
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Among the singularities of the linear ODE for <J>^ given in (|8.2[) . (|8.3p or obtained from the 
formula (29) given in [30J up to n = 15, we have found no other singularity identified with the 
remarkable Heegner numbers [100] or, more generally, with other selected values of the modular 
j-function, associated to complex multiplication. 

Could it be that the (non-Nickelian) singularities (|8.2p . (|8.3p . which do not match with 
complex multiplication of the elliptic curves, are actually selected for mathematical structures 
more complex or more general than elliptic curves (possibly linked [39] to SL(3, Z) instead of 
SL(2, Z) modular group)? This could amount to moving away from the isotropic Ising model 
towards the Baxter model. At first sight the analysis of the anisotropic Ising model [26] could be 
considered as a first step in that "Baxter-model" direction. The selected situations for elliptic 
functions and complete elliptic integrals, would thus, be generalized to the search of "selected 
situations" of their multidimensional generalizations (Lauricella, Kampe de Ferie, Appell, . . . ) 
that we have actually seen to occur in the anisotropic Ising model |26j and even in our series 
expansions of x^ an d X^ ■ 

Along similar lines, one may recall the n-fold integrals introduced by Beukers, Vasilyev [1161 
[TT7] and Sorokin pm EDI 



i a +l (1 - tj)Z 

dtidt 2 ■ • ■ dt a+ i | 

r- 



(z - ht 2 ■ ■■ta+i) 3 
and other well-poised hypergeometric functions 



((2r + l)n + 2)! fr+1 



n |2r+l 
/ 



(r+l)n+l 



JO 

n 



a+1 \ 
7=1 



(z-t 1 t 2 ---t a+1 ) 2r + 1 



V 



dtidt 2 ■ ■ ■ dt a+ i 



z + t\t 2 ■ ■ ■ t a+ i 

Z - tit 2 ■ ■■t a+ l) 3 ' 



J 



or the Goncharov-Manin integrals [13] which occur in the moduli space of curves |19l US] . These 
integrals [311 HH H21 [551 EH EH1 1101 j look almost the same as the ones we have introduced and 
analyzed in the study of the diagonal susceptibility of the Ising model [5B] . 

It is worthy to recall that ("(3) appeared in some of our "connection matrix method" results 
for the differential Galois group |128| of the Fuchsian linear ODE for an d X^\ an d the 
occurrence of zeta functions in many n-fold integrals. Also recall that Feynman amplitudes can 
be seen as periods in the "motivic sense" [19], and are often linked to multiple zeta numbers. 
Along this line, the following integral [411 143| deals with C(3): 



In{z) 



1 r (1 -u) n u n (l-v) n v n (l-w) n w n 

dudvdw - 



(1 



uv 



in+1 



uvw 



in+1 



(9.13) 



From the series expansion of this holonomic n-fold integral, we have obtained [30] an order 
four Fuchsian linear differential equation (see Appendix[F]). On such linear differential operators 
the "logarithmic" nature of these integrals becomes clear. The occurrence of linear differential 
operators is not a complete surprise if one recalls that in Apery's proof of the irrationality of £(3) 
a crucial role is played by the linear differential operator [18] 



L = (t 2 - 34t + l)t 2 £> f 3 + (6t 2 - 153t + 3)tD? + {It 2 - 112* + l)D t + (t - 5) 
this operator being linked to the modularity of the algebraic variety: 



x H h y + 

x y 



1 1 1 

- + z + - + w + — = 0. 

V z w 
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These n-fold integrals have to be compared with the (more involved) re-fold Ising integrals 
corresponding to the x^ n \ an d to the theory of elliptic curves (rather than rational curves CP\ 
in the previously cited examples (3TJ [3U HH H21 [551 ESI [691 HOI] ) , we try to underline in this 
paper. 

With these new singularities, are we exploring some remarkable "selected situations" of some 
moduli space of curves [H [32] corresponding to pointed [91 [38l EI] (marked) curves [10] , instead 
of simple elliptic curves [53] ? In practice this will probably just correspond to considering a pro- 
duct of n times a rational or elliptic curve minus some sets of remarkable algebraic varieties [28] , 
X{Xj = 1, XiXjXk = 1, hyperplanes Xi = Xj, .... 

10 Conclusion 

We have displayed several examples of n-fold holonomic integrals associated with the two- 
dimensional Ising model on a square lattice [8]. The corresponding linear differential opera- 
tors with polynomial coefficients are shown to be very closely linked to the theory of elliptic 
curves (and modular forms) and display many remarkable structures (Russian-doll structure, 
direct sum structure, complex multiplication as selected singular values for these operators, . . . ). 
These linear differential operators are not only Fuchsian operators, they are Fuchsian operators 
with rational exponents: the various indicial polynomials corresponding to all the regular singu- 
larities of these linear differential operators have only rational (or integer) roots. It is tempting 
to try to understand these deep algebraico-differential structures as a consequence of the under- 
lying elliptic curve in the Ising model, or more generally, of some algebraic varieties built from 
this elliptic curve (product of curves, . . . ) , or corresponding to the integrands of these n-fold 
integrals. Could it be possible that these large number of remarkable properties have a geo- 
metrical interpretation (generalisation of hyper geometric functions and Picard-Fuchs systems, 
Griffiths cohomology of hypersurface of CP n , rigid local systems [171 SSI E21 (S3 EH ESI 166] [99] . 
. . . ) with a strong background of algebraic geometry? One could, for instance, imagine that 
these various n-fold holonomic integrals might be interpreted as periods of some algebraic vari- 
eties, all the strong and deep algebraico-differential structures we have displayed in this paper, 
being a consequence of this very rigid geometrical framework. The central role played by the 
theory of elliptic curves and their isomonodromic deformations (Painleve equations) for the Ising 
model on a lattice is also underlined in the fundamental finite-difference (non-linear overdeter- 
mined) system of quadratic functional relations |82[ \M\ [96"1 [97] (see (jA.ip in Appendix [A]) for 
the two-point correlation functions of the Ising model on the square lattice. As Painleve and 
(discrete) integr ability specialists call it, these lattice equations are finite-difference generalisa- 
tion of Painleve equations and they have a lot of very deep consequences : they are, for instance, 
the very reason why the susceptibility series can be calculated from a program with polynomial 
growth [93] . Such an overdetermined system (jA.ip can be seen as generating an infinite number 
of non-trivial identities on the complete elliptic integrals of the first and second kind. 

It is important to note that all these remarkable structures and deep symmetries (remarkable 
functional identities, algebraico-differential structures, modular forms, continuous [811 [90] and 
discrete Painleve structures, . . . ), underline the central role played by the theory of elliptic curves 
for the two-dimensional Ising model on a lattice. Note that a large part of these remarkable struc- 
tures and deep (lattice) symmetries is lost in the scaling limit. In the scaling limit some of these 
remarkable structures remain (the Russian-doll telescopic embedding of the linear differential 
operators), but, for instance, the direct-sum structure is lost. The scaling limit yields the occur- 
rence of an irregular singularity at infinity: the Fuchsian character of the linear operators is lost, 
as well as most of the remarkable structures associated with the underlying elliptic curve theory. 
For instance for two-point correlation functions, the complete elliptic integrals of first and sec- 
ond kind K and E are replaced by modified Bessel functions (with their irregular singularity at 
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infinity), but the fact that form factors are simple polynomial expressions of E and K is lost: 
the form factors, in the scaling limit, are not simple polynomial expressions of modified Bessel 
functions. In the scaling limit, a large part of the strong background of algebraic geometry that 
exists on the lattice model, and yields so many remarkable deep and strong structures and sym- 
metries, seems to disappear. If the geometrical interpretation we suggested for the lattice model 
exist, could it be possible that it is essentially lost in the scaling limit, the underlying algebraic 
varieties necessary for this geometrical interpretation being lost, or becoming some complicated 
analytical manifolds? Recalling the emergence of an irregular singularity (at infinity), an irregu- 
lar singularity can, in principle, be understood [2] as a confluence of two regular singularities 
(for complete elliptic integrals of first and second kind we have the confluence of the two regular 
singularities 0, oo among the three regular singularities 0, 1, oo). To our knowledge we have not 
often seenH in the litterature the structures associated to irregular singularities (Stokes multipli- 
ers, singular behaviours, . . . ) be obtained as a "confluent limit" of the structures associated with 
the two regular singularities. From a general viewpoint, in a desire to see analytical manifolds as 
a confluent limit of algebraic varieties, one can imagine that the structures of the Ising model in 
the scaling limit could, in principle, be obtained from a (very involved) "confluent limit" of the 
remarkable structures deeply linked to the theory of elliptic curves that exist for the Ising model 
specifically on the lattice. This remains to be done. In practice we see that, paradoxically from 
a criticality-universality mainstream viewpoint, the (off-critical, non-universal) Ising model on 
a lattice has much deeper, and fundamental, structures than the Ising model in the scaling limit. 
Note that the results of holonomic and algebraic-geometry nature we have displayed in this 
paper, are not specific of the two-dimensional Ising model, or, even, of free-fermion models. We 
have not used the free-fermion character of the Ising model. We have heavily used the elliptic 
parametrisation of the two-dimensional lattice Ising model. One can imagine that many of these 
results, and structures, exist for Yang-Baxter integrable models with an elliptic parametrisa- 
tion (the Baxter model \12\ [23l I24j . . . . ), and, more generally, for any Yang-Baxter integrable 
modej^l, the central role of the elliptic curve being replaced by the relevance of the algebraic 
variety emerging in the Yang-Baxter equations (higher genus curves [5], Abelian varieties, . . . ). 
Integrable models on a lattice are probably deeper, and dressed with much more symmetries 
and remarkable structures^!, than their scaling limits. Such an apparently paradoxical (for the 
field theory mainstream) conclusion is certainly not a surprise for Painleve and (discrete) inte- 
grability specialists who are used to see, and understand, lattice equations as deeper, and more 
fundamental [85 j 1106] . than the differential equations. 

A Quadratic partial difference Painleve generalisations 

Quadratic partial difference equations were shown [82j [MJ ISSl EZ] to be satisfied by two-point 
correlation functions of the two-dimensional Ising model on the square lattice. These quadratic 
partial difference equations (valid in the anisotropic case), are actually valid for the A-extension 
of the two-point correlation functions C(M, N; A) for any value of X : 

s 2 [C{M, N; A) 2 - C(M, N - 1; A)C(M, N + l; A)] 

+ [C*{M,N;X) 2 -C*(M - 1, N; X)C* {M + 1, N; X)} =0, (A.l) 
s' 2 [C(M, N; A) 2 - C(M — 1, N; X)C(M + 1, N; A)] 

23 With the exception of C. Zhang explaining |130] Ramis' confluent approach [471 11021 1112| of irregular singu- 
larities. 

24 It is not even clear that one has to restrict oneself to Yang-Baxter integrability: n-fold integrals associated 
with particle physics, polyzeta functions, Feynman diagrams, etc . . . seem to indicate that one has a much more 
general framework for these ideas [T01 172] . 

25 See for instance Baxter's concept of Z-invariance [2l 112], 
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+ [C*{M, N; A) 2 - C*{M, N - 1; X)C*(M, N + 1; A)] = 0, 
ss'[C(M, N; A)C(M + 1, N + 1; A) - C(M, JV + 1; A)C(M + 1, N; A)] 

= C*(M, N; X)C*(M + 1, JV + 1; A) — C*(M, N + 1; A)C*(M + 1, N; A)] 



with s = sinh(2K), s = sinh(2i<r), where and i^' correspond to the horizontal and vertical 
coupling constant of the anisotropic square Ising model and where C* are the A-extension of the 
dual correlation functions. 

B Algebraic solutions of PVI for A = cos(7rra/n) 
and modular curves 

The unexpectedly simple expressions for the form factors f^\ of Sections [2HU and the corre- 
sponding remarkable differential structures, may be used to obtain many further results. We 
displayed some of these results in Section 6 of [25]. Recalling that, when A = 1, the Ising 
correlation functions C(N,N; 1) satisfy Fuchsian linear differential equations [27] with an order 
that grows with N, it is quite natural to inquire whether there are any other values of A for 
which C(N,N;X) will satisfy a Fuchsian linear differential equation. One such family of A is 
motivated by the work of Cecotti and Vafa [33] on N = 2 supersymmetric field theories where 
they encountered A extensions of the Ising correlations in the scaling limit [83] with (m and n 
are integers): 

A = cos(7rm/n). 

Indeed, we have found that for n = 3, . . . , 20, the functions C(N, N; A) satisfy Fuchsian linear 
differential equations whose orders, in contrast with those of the A = 1 equations [27], do not 
depend on N. 

The function C(N,N;X) is such that its log-derivative is actually a solution of the sigma 

form of Painleve VI : it is a transcendent function "par excellence" . However, the unexpectedly 

(i) 

simple expressions for these form factors /^ jv» strongly suggest to try to resum the infinite 

sums (11. 5ft . and (II. 6ft . of form factors fjp N , corresponding to the function C(N, N; A), and see 
if these transcendent functions could be "less complex" than one can imagine at first sight, at 
least for a set of "singled-out" values of A. For instance, are there any values of A ^ 1 which 
share, with A = 1, the property that C(N, N; A) satisfies a Fuchsian linear differential equation? 

Actually, introducing, instead of the modulus k of elliptic functions (for T > T c , k = s 2 ), or 
the s and t variables, the nome of the elliptic functions (see relations (5.7)-(5.11) in |93j), we have 
been able to perform such a resummation, getting, for arbitrary A, nice closed expressions for the 
C(N, N; A) for the first values of N, (N = 0, 1, 2, . . . ), as sums of ratios of theta functions (and 
their derivatives), corresponding to Eisenstein series, or quasi-modular forms. These results will 
be displayed in forthcoming publications. The simplest example corresponds to = where 
C-(N, N; A) is just the ratio of two Jacobi 83 functions 



All these results strongly suggest to focus on u = nm/n (m and n integers) yielding for the 
possible choice of "selected" values of A: 



C_(0,0;A) 



03(0,9) 



where 



A = cos(u). 



(B.l) 



A = cos(7rm/n). 



(B.2) 



Actually these special values (|B.2p of A already occurred in a study of = 2 supersymmetric 
field theories [33] in a similar series construction of solutions of the Painleve V (or Painleve III 
for a ratio of functions) equation for the scaling limit of the Ising model [123] . 
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Recalling the quadratic finite difference equations [3 [M] (| A. 1|) we can deduce that the off- 
diagonal terms C(M,N;X) are, in the isotropic case, algebraic expressions of sum of ratios of 
theta functions and their derivatives. For the singled-out values A = cos(irm/n) (n and m 
integers), the off-diagonal terms C(M,N;X) are, in the isotropic case, algebraic expressions of 
the variable t: do these algebraic expressions also correspond to modular curves? Actually they 
clearly single out t = 0,1, oo, .... 

Remark. For this set of selected values of A the A-extension C(N, N; A) are seen to be algebraic 
expressions of the variable t and, more remarkably, associated with a modular curve P(C, t) = 
(P denotes a polynomial with integer coefficients, C denotes C(N, N; A) for A = cos(7rm/n), for 
certain integer values of n and m, and the only branching points are for t = 0, 1, oo). The fact 
that the only singular points are t = 0, 1, oo can be seen to be inherited from the fact that the 
A-extension C(N,N;X) is actually solution of (|1.7j) for any X: the sigma form of Painleve VI, 
namely (11 . T[> . naturally singles out t = 0, 1, oo (and only t = 0, 1, oo). 

C Linear differential equations of some <f>^ 
C.l Linear ODE for 

(3) 

The minimal order linear differential equation satisfied by & H reads 

n=0 

where 

a 5 (w) = (l-w)(l- Awf (1 + Awf (1 + 2w) (1 + 3w + iw 2 )w 3 P 5 (w), 
a 4 (w) = (1 - 4w) 3 (1 + Aw) w 2 P 4 {w), a 3 {w) = -2 (1 - Aw) 2 wP 3 (w), 
a 2 (w) =2(1- Aw)P 2 {w), ai {w) = -8Pi(u>), a (w) = -96P (w), 

and where the apparent singularities polynomial P${w) in the head polynomial reads 

P h (w) = -5 + 2lw + 428w 2 + 5364w 3 - 82416w 4 - 299504w 5 + 714944™ 6 

+ 3127872w 7 - 8220672-u; 8 - 25858048w; 9 - 7077888u; 10 + 31424512m; 11 
- 42467328-u; 12 - 31457280w; 13 - 4194304u; 14 + 4194304u; 15 , 

the other polynomials P n (w) are given in [30]. 



C.2 Linear ODE for 



The minimal order linear differential equation satisfied by & H reads (with x = 16w 2 

6 d n 
^>" (X) ^ F(X)=0 ' 

n=0 

where 

a 6 (x) = 64 (x - 4) (1 - x) A x A P 6 (x), a 5 (x) = -128(1 - x) 3 x 3 P 5 (x), 

a A (x) = 16(1 - x) 2 x 2 P A {x), a 3 (x) = -64(1 - x)xP 3 (x), 

a 2 {x) = -AP 2 {x), ai {x) = -8Pi(x), a {x) = -3(1 - x)P (x), 
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where the apparent singularities polynomial Pq(x) in the head polynomial reads: 

P 6 (x) = 128 + 2233x - 2847x 2 + 3143rr 3 - 3601x 4 + 144x 5 - 64x 6 . 
The other polynomials P n {w) are given in [30] . 

C.3 Linear ODE modulo a prime for <& H 

(5) 

The linear differential equation of minimal order seventeen satisfied by <& H is of the form 

n=0 

with 

a 17 (w) = (1 - Aw) 12 (1 + Aw) 9 (1 - wf (w + 1) (1 + 2w) (l + 3w + Aw 2 ) 2 

x (1 - 3w + w 2 ) (1 + 210- Aw 2 ) (1 + Aw + 8tt; 2 ) (l — 7w + 5w 2 - 4w 3 ) 
x (1 - to - 3w 2 + Aw 3 ) (l + 8w + 20w 2 + 15w 3 + iw A ) w l2 P 17 (w), 

a 16 (w) = w 11 (1 - 4w) u (1 + iw) 8 (1 - w) (l + 3w + 4w 2 ) PieH, 

ai 5 («;) = u; 10 (1 - 4w) 10 (1 + 4u;) 7 P 15 (w), 

a u (w) = w 9 (l- Aw) 9 (1 + Awf Pu(w), 



where the 430 roots of Pn(w) are apparent singularities. The degrees of these polynomials P n (w) 
are such that the degrees of ai{w) are decreasing as: deg(oj + i(w)) = deg(aj(w)) + 1. In fact, 

((5) 

with 2208 terms we have found the ODE of ' at order q = 28 using the following ansatz for 
the linear ODE search {Dw denotes d/dw) 

z(i)p(i)Dw l 

i=0 

with 

z (i) = w a (~ 1+i \l - I6w 2 ) a ^ 1+i h« {1+t ~ q) 

where a(n) = min(0, n) and: 

zq = (1 + w)(1 -w)(1 + 2w)(1 -3w + w 2 ){\ + 2w- Aw 2 ){\ + 3w + Aw 2 ) 
x (1 + Aw + 8w 2 )(l -7w + 5w 2 - Aw 3 )(l - w - 3w 2 + Aw 3 ) 
x (1 + 8w + 20-u; 2 + 15w 3 + Aw 4 ) 

the p(i) being the unknown polynomials. The minimal order ODE is deduced from the set of 
linearly independent ODEs found at order 28. Instead of these linear ODEs with quite large 
apparent singularities polynomials, we can provide an alternative linear ODE of higher order 
with no apparent singularities. This is the so-called "desingularization" procedure of a linear 
ODE. The price to pay to get rid of the large apparent polynomial can be that the higher 
order ODE with no apparent polynomial may not be Fuchsian anymore (because of an irregular 
singularity at infinity). One can also consider desingularizations preserving Fuchsianity. 
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We give below the linear differential operator of order 25, modulo the prime 27449. At this 
order, the linear differential operator has no apparent singularities 

(w + 12523) 2 (™ + 3989) (to + 9241) (w + 21789) 2 (w + 20587) 9 
x (™ 3 + 6861™ 2 + 6864™ + 6862) (w + 6862) 12 (w + 27448) 2 
x (™ 3 + 14565™ 2 + 6775™ + 7627) (w + 13725) 
x (™ 3 + 20586™ 2 + 6862™ + 20587) (w + 4483) (w + 8965) 
x (w + 19750) (w + 18481) (w + 1) (w + 9736) w 12 Dw 25 
+ 1134 (w + 12523) P23 (w + 21789) (w + 20587) 8 
x (™ + 6862) 11 (w + 27448) w u Dw 24 

+ 24190 (w 2 + 26576™ + 1920) (w + 6798) P 17 (w + 20587) 7 

x (™ + 6862) 10 P 10 P9<59^ 10 ^ 23 

+ I8886P7 (to + 16985) (w + 20587) 6 (to + 6862) 9 P 22 

x (™ 4 + 1460™ 3 + 5534™ 2 + 16322™ + 22102) Q 7 P n P 51 w 9 Dw 22 

+ 2514™ 8 P 82 (w 2 + 24026™ + 4328) (w + 20587) 5 P u 

x (™ 6 + 834™ 5 + 1985™ 4 + 9290™ 3 + 16354™ 2 + 4746™ + 4283) 

x (to + 6862) 8 (w + 10153) Dw 21 

+ 7703P 24 ™ 7 (w 2 + 21361^ + 14160) (™ + 20587) 4 P 80 
x (to + 6862) 7 (w + 12463) Dw 20 

+ 7921Pi4u; 6 (w 2 + 9437™ + 10364) P 7 (w + 20587) 3 P 62 
x (™ + 4208) Pio (w + 3899) Dw 19 

+ 26675 (w 2 + 7849™ + 652) P 61 (to + 20587) 2 (u> + 6862) 5 P 42 

x (™ 6 + 16220™ 5 + 14100™ 4 + 16063™ 3 + 18759™ 2 + 27385™ + 22000) 

x (™ + 6862) 6 Pi 2 ™ 5 L>™ 18 

+ 6905™ 4 P 9 (w + 20587) (w + 6862) 4 Q 9 P 71 (w + 21055) (w + 5262) 

x (™ 5 + 24655™ 4 + 5244™ 3 + 12152™ 2 + 16121™ + 18032) P 17 Dw 17 

+ 12508 (w + 10977) w 3 (w + 6862) 3 (™ 3 + 8559™ 2 + 11289™ + 9639) 

x (™ 3 + 14200™ 2 + 21840™ + 15970) 

x (™ 5 + 2056™ 4 + 12162™ 3 + 6126™ 2 + 22729™ + 16548) 

x P16P77 (w 2 + 16577™ + 6876) P 8 Dw lG 

+ 8317™ 2 (™ + 21290) (™ 2 + 4503™ + 23586) 

x (™ 5 + 26400™ 4 + 26455™ 3 + 17876™ 2 + 7063™ + 15487) 

x (™ 6 + 10703™ 5 + 1763™ 4 + 12902™ 3 + 6843™ 2 + 2123™ + 18674) 

x P43P45 (w + 17888) (w + 6862) 2 Pi 3 Dw 15 

+ 13440P 68 ™ (™ + 6862) P 48 (w + 13354) £>™ 14 + 11221P 27 P 39 P 8 Q 22 P 22 L>™ 13 
+ 4954 (™ + 15012) P 46 (™ 2 + 15293™ + 14779) P 68 £>™ 12 
+ 26649Pi 06 P8 (w + 11457) (w + 24683) Dw 11 

+ 20358 (™ 2 + 2503™ + 22492) P 22 (™ 3 + 13156™ 2 + 2684™ + 8942) P 7 
x (™ 4 + 15516™ 3 + 4945™ 2 + 4117™ + 4952) P 77 Dw w 
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+ 18965 (w 2 + 10905u; + 26489) 

x (w 5 + 13584w; 4 + 17617-u; 3 + 21809w 2 + 8787m; + 16558) PiiP A0 P^D 
+ 16238 (w 5 + 6787u; 4 + 23738w 3 + 20731w 2 + 22246^ + 21284) 
x (w 6 + 13024w 5 + 11465w 4 + 2372w 3 + 6908w 2 + 612w + 10791) 
x (w + 12744) P wl Dw 8 + 706 (w + 3771) P^P A& Dw 7 
+ 16695 (w 5 + 8888w 4 + 6808w 3 + 26162u; 2 + 12699w + 21123) 
x (w 3 + 14297w 2 + 5276w + 11625) P 79 (w 2 + 17252u; + 11894) 
x (w 5 + 1644w 4 + 15551w 3 + 18064w 2 + 8721w + 21521) P l7 Dw 6 
+ 19013P 2 4^3i (w 3 + 4946u> 2 + 25614w + 25774) 
x (w 5 + 884u; 4 + 16517u; 3 + 2210u> 2 + 13068u; + 15270) P 37 P w Dw 5 
+ 12971 (w 2 + 7070w + 2438) P 25 

x (w 6 + 869w 5 + 26611u; 4 + 21387w 3 + 25750w 2 + 21941m + 24186) 
x (w + 23774) P 17 P m P A2 Dw 4 

+ 19579P 2 9^n (w + 19783) (w + 9767) P 7 P 36 (w + 15147) P 22 Dw 3 

+ 9934 (w 6 + 7835u) 5 + 21707m 4 + 5095m 3 + 8057m 2 + 13773m> + 14223) 

x (w 4 + 8474m 3 + 901 lw 2 + 12522m; + 22902) P w P 72 P 9 Dw 2 

+ 1875Q 2A P 51 P 7 P 2A Dw + 4564 (w + 25617) P 104 , 

where the P n and Q n is a "short" notations to encounter polynomials of degree n (they can be 
different from one coefficient of the m-th derivative Dw m to another). 

Remark. We sketch such a quite tedious result (if we give explicitly the undefined polynomials 
the result would be really huge . . . ) to give the reader some hint of how such an exact result 
modulo a prime looks like: the exact expressions of the various polynomials, which are the 
coefficients in front of the derivatives, can actually be factorized modulo prime without any 
ambiguity. For instance the factor (w + 27448) 2 in the head polynomial (coefficient of Dw 25 ) is 
nothing but (w — l) 2 modulo the prime 27449. The interest of such an exact calculation is that 
we can exactly compare the various factors in the head polynomial with a set of polynomials we 
have conjectured to be singularities of the linear ODE. We can totally confirm the existence of 
some of (or all) these conjectured polynomials, and discriminate between apparent singularities 
and "true" singularities. The prime is large enough to avoid any ambiguity corresponding to 
accidental factorisations (because the prime would be too small): this is confirmed by the same 
calculations performed for other similar large enough primes. 

C.4 Linear ODE modulo a prime for & H 

The linear differential equation of minimal order (namely twenty-seven), satisfied by $^ , reads 
(with x = w 2 ) 

27 d n 

n=0 

with 

a 27 { x ) = (1 - 16x) 16 (1 - Axf (1 - x) (1 - 25x) (1 - 9x) x 21 
x (1 - x + 16x 2 ) (1 - lOx + 29x 2 ) P 27 (x), 
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026 0) 


= 0- 


16x) 15 


(l-4x) 2 x 20 P 26 O), 


025 0) 


= 0- 


16x) 14 


(l-4x)x 19 P 25 (x), 


£124 0) 


= 0- 


16x) 13 


X 18 P 2 40), 



where the 307 roots of -P27O) are apparent singularities. The degrees of the P n (w) polynomials 
are such that the degrees of ai(w) are decreasing as: deg(aj+i(w;)) = deg(oj(w;)) + 1 

In fact, with 1838 terms we have found the linear ODE of & H at order q = 42 using the 
following ansatz for the linear ODE search {Dx denotes d/dx) 

z(i)p{i)Dx l 

i=0 

with 

z (i) = x<~ l+l \l - 16x) a ^ 1+l h^ 1+l - q) , 

where a(n) = min(0, n) and 

zq = (1- 25s)(l - 9x)(l - 4x)(l - x)(l -x + 16x 2 )(l - Wx + 29x 2 ), 

the p{i) being the unknown polynomials. The minimal order ODE is deduced from the set of 
linearly independent ODEs found at order 42. 

Here also, instead of this linear ODE with quite a large apparent singularities polynomial, 
we can provide an alternative linear ODE of higher order with no apparent singularities (but 
it may not be Fuchsian anymore). We give in the following the linear differential operator, 
modulo the prime 32749, of order 30. At this order, the linear differential operator has no 
apparent singularities 

(x + 8187) 3 2 + 10234x + 22515) (x + 10234) 16 (x + 32748) x 21 
x + 10113) + 31439) + 14555) (x + 15860) Dx 30 
+ 8594 (x + 8187) 2 (x + 10234) 15 (x + 551) x 20 
x (x 4 + 22258x 3 + 24734x 2 + 32441x + 31408) P 85 P 2 iDx 29 
+ 30840 (x + 8187) (x + 10234) 14 (x + 28552) P b0 P & i 
x (x 3 + 5377x 2 + 7622x + 28946) x l9 Dx 28 
+ 22635 (x + 10234) 13 x 18 (x + 22434) P 16 P 23 P m P 8 Dx 27 
+ 15369 (x + 10234) 12 (x + 3) (x + 25968) (x + 15827) x 17 
x (x 6 + 704x 5 + 19667x 4 + 16573x 3 + 221x 2 + 5237x + 24649) 
x P 7 P 8 P 43 P 2 2PwPi3Dx 26 

+ 3485 (x + 10234) 11 Pi 09 x 16 (x 5 + 6389x 4 + 9765x 3 + 14807x 2 + 31264x + 20696) 
x (x 4 + 24125x 3 + 31305x 2 + 16748x + 12080) (x + 30475) Dx 25 

10 15 I J2 1 /hqo„ 1 td T3 D__ r>™24 



+ 30663 (x + 10234) 10 x 15 (x 2 + 4182x + 14901) P 9 P A9 P m D 



x 



100 



+ 286 (x + 10234) 9 x 14 (x 3 + 15593x 2 + 11835x + 28482) P 
x (x 3 + 3549x 2 + 28115x + 25784) (x 3 + 26209x 2 + 12548x + 5267) P 12 Dx 23 
+ 3836 (x + 10234) 8 P 28 P 29 P m P 7 (x + 495) x 13 £>x 22 

+ 16272 (x + 10234) 7 (x + 22272) x 12 (x 2 + 10657x + 16936) P 28 P 7& Pi 6 Dx 2 



.21 
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+ 108 (x + 10234) b x 11 (x + 12405) (x + 17141) (x + 16254) P u P 47 Dx M 
+ 29777 (x + 10234) 5 x 10 P 23 P 15 P n Pm 

x (x 6 + 23995x 5 + 4070x 4 + 27561x 3 + 19739x 2 + 3632x + 18638) 

x (x 2 + 24371x + 11409) (x + 22618) 

x (x 4 + 6026x 3 + 10330x 2 + 10566x + 27129) 

x (x 3 + 1171x 2 + 10654x + 2741) Px 19 

+ 6325 (x + 10234) 4 (x + 20731) x 9 P 42 P25P51 {x 2 + 16210x + 27515) 
x (x 5 + 11144x 4 + 6536x 3 + 25134x 2 + 10963x + 29010) Dx 18 
+ 11986 (x + 10234) 3 (x 2 + 16634x + 15614) P m P 13 P 72 x s Dx 17 
+ 16154 (x + 10234) 2 (x 3 + 28948x 2 + 8126x + 18460) x 7 
x (x + 12041) (x + 9774) P 4S PuP3iDx w 
+ 5724 (x + 10234) x 6 P 25 P 99 

x (x 4 + 5795x 3 + 12069x 2 + 26629x + 1320) (x + 17629) Dx 15 
+ 9093^10^39^34^45 {x 2 + 24077x + 23664) x 5 Px 14 
+ 7105x 4 Pio 5 P 2 2 {x 3 + 1115x 2 + 1326x + 29632) Dx 13 

+ 15661 (x 5 + 25773x 4 + 9200x 3 + 26470x 2 + 25643x + 1121) P 117 P 8 x 3 Dx 12 

+ 15107 (x 4 + 14631x 3 + 6554x 2 + 7715x + 3048) (x + 6577) 

x (x 3 + 31616x 2 + 26256x + 9612) x 2 

x (x 3 + 4756x 2 + 28396x + 28874) P l0 P b %P bl Dx 11 

+ 26871 (x 2 + 11848x + 20401) PnP m P m xDx w 

+ 541P 8 g 8 Pi4P66^33 {x + 932) Dx 9 + 4081 (x 2 + 6617x + 3717) PisPmPisPm 

x (x 3 + 671x 2 + 11514a; + 23683) (x 2 + 18485x + 5460) 

x (x 6 + 23400x 5 + 12243x 4 + 21913x 3 + 27012x 2 + 17751x + 12915) 

x (x 2 + 21416x + 31601) P 7 Px 8 + ^lP 9 P m PuP 71 P 7 Dx 7 

+ 6649P 73 P 2 6P23 (x 3 + 4619x 2 + 29249x + 15768) (x 2 + 12299x + 30824) Px 6 

+ 2785 (x + 16825) (x + 27878) (x + 15523) P7P9P102 

x (x 5 + 13600x 4 + 24394x 3 + 31753x 2 + 19488x + 21782) Dx 5 

+ 14219Pi 8 Pio5 (x 2 + 19404x + 23792) Px 4 

+ 28992Pn7 (x 3 + 9172x 2 + 23091x + 20852) (x + 15969) 

x (x 2 + 31441x + 5617) (x + 32359) Dx 3 

+ 23799P 7 Pio2 {x 3 + 7169x 2 + 15714x + 24623) 

x (x 6 + 21946x 5 + 31546x 4 + 16591x 3 + 19174x 2 + 23949x + 23818) 

x (x + 10738) (x 4 + 10798x 3 + 20551x 2 + 6303x + 7193) Dx 2 

+ 3515 (x + 16461) (x + 6204) Pi 5 Pi P 8 4 (a; + 4637) (x + 1197) 

x (x 6 + 31051x 5 + 11003x 4 + 821 lx 3 + 28599x 2 + 20034x + 25604) 

x (x + 17071) (x 2 + 26670x + 1134) Px + 25380 (x + 30713) P35P15P33P26 

x (x 5 + 8267x 4 + 15086x 3 + 11158x 2 + 26216x + 31098) 

x (x 6 + 30193x 5 + 28390x 4 + 17930x 3 + 26696x 2 + 7578x + 16219) , 
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where the P n and Q n a "short" notation for polynomials of degree n (that may be different from 
one order D" 1 to another). The factor (x + 32748) in the head polynomial (coefficient of Dx 30 ) 
is nothing but the factor (x — 1) modulo the prime 32749. 

(7) (8) 

D Singularities in the linear ODE for $> H and $> H 

(7) 

For <J>^ y , we generated long series, unfortunately, insufficient to obtain the corresponding linear 
ODE. Actually, we have also generated very long series modulo a prime (40000 coefficients) 
and we have not been able to find a linear ODE when the order of the ODE is less than 100. 
However, by steadily increasing the order q of the ODE and the degrees n of the polynomials in 
front of the derivatives, one may recognize, in floating point form, the singularities of the linear 
ODE as the roots of the polynomial in front of the higher derivative. A root is considered as 
singularity of the still unknown linear ODE, when as q and n increase, it persists with more 
stabilized digits. 

(7) 

Using 1250 terms in the series for & H , the following singularities are recognized 

(1 - Aw) (l - 5w + 6w 2 - w 3 ) (l + 2w - 8w 2 - 8w 3 ) (1 + Aw) (l + w)w 

x (1 + 2w - w 2 - w 3 ) (1 - 3w + w 2 ) (1 + 2w - Aw 2 ) (l - 7w + 5w 2 - Aw 3 ) 
x (l-3w- 10w 2 + 35w 3 + 5w 4 - 62w 5 + 17w; 6 + 32w 7 - 16w 8 ) 
x (1 + 8w + I5w 2 - 2lw 3 - 60w 4 + 16u> 5 + 96w e + 6Aw 7 ) 

x (l - Aw - 16w 2 - A8w 3 + 32u; 4 - 128w 5 ) (l - lOw + 35w 2 - 51w 3 + 2lw A - Aw 5 ) 
x (1 + 7w + 26w 2 + 7w 3 + Aw 4 ) (l + 8w + 20w 2 + 15w 3 + Aw 4 ) 
x (1 + 12™ + 5Aw 2 + I12w 3 + 105w 4 + 35w 5 + Aw 6 ) = 0. 

Note that we have not seen with the precision of these calculations the occurrence of the 

(3) 

singularities of the <& H . 

(8) 

With similar calculations using 2000 terms for $Vr , the following singularities are recognized 

(1 - 2w) (1 + 2w) (1 - 2w 2 ) (1 - Aw) (l - Aw + 2w 2 ) (1 + Aw) (1 + 3w) w 

x (1 + Aw + 2w 2 ) (1 - 8w 2 ) (1 - 3w) (1 - w) (1 + w) (1 - 5w) (l + 2w 2 ) 
x (1 - 26w 2 + 2A2w 4 - 960w 6 + 1685u; 8 - 1138u; 10 ) (l - 10w 2 + 32-u; 4 ) 
x (1 - 30w 2 + 56w 4 - 1312w 6 ) (l - 6w + 10w 2 ) (1 - 6w + 8w 2 - Aw 3 ) 
x (1 + 5w) (l + 6w + Ww 2 ) (l + 6w + 8w 2 + Aw 3 ) = 0. 

Note that the stabilized digits in these singularities can be as low as two digits. 



E Selected values for Liouville theory and Potts models 

New classes of critical statistical models where suggested [46] by Gervais and Neveu from the 
construction of Liouville field theory. With the Q-state standard scalar Potts model notations 
(see (1.3) in [46]), they introduced y, such that Q 1 ! 2 = 2cos(7ry/2). Rational values of y 
correspond to selected values of Q (Tutte-Beraha numbers see Section 4 of |104j ) for which 
the standard scalar Potts model has rational critical exponents. At this step, and in order to 
make explicit the selected role of these particular values, we can recall the expression (see (3.3) 
in [103] ) of the partition function per site of the Q-state standard scalar Potts model on the 
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checkerboard lattice in terms of Eulerian products (see (3.5) in [103] ) like (with the notations 
of [103] ): 



= n t^Ws' where Ql/2=t+ i (K1) 

n=l 

This Eulerian product form made very clear the fact that the partition function can be seen as 
some automorphic function with respect to an infinite discrete group generated by the inverse 
relation and the symmetries of square [57[ [58] . Such Eulerian product over an infinite discrete 
group also made very clear the fact that these singled-out values of Q actually correspond 2 "^! to 
iV-th root of unity situation 

t N = 1, where Q 1/2 = t + j (E.2) 

that occur in some many domains of theoretical physics |111[ I113| (dilogarithms, Kac determi- 
nant, . . . ). Do note that such situation generalizes, mutatis mutandis, to the Baxter model: the 
partition function per site can actually be written as an infinite discrete product pTJ[13j[T3] over 
a group generated by the inverse relation and geometrical symmetries of lattice [76] , expressions 
like (jE.ip being replaced by (with Baxter's notations [TTJ [131 1 

00 00 - „m„.4n-l. 



= n n f 

m=0 n=l 



-.mj.4n-\-lg ' 

m=U n=l 

where 

q = exp(— 7rX'/T), x = exp(— ir\/2Z), z = exp(— nv/2I). 

Such an expression of the partition function per site of the Baxter model as infinite product 
can also be found in [39J in terms of product and ratio of theta and elliptic gamma functions. 

In [36] Gervais and Neveu underlined that they had built Liouville field theory for other 
singled-out values of Q than iV-th root of unity situations like (|E.2j) . namely (see (2.3) in 



Q = -4 sinh 2 (71-/2^3), and Q = 4cosh 2 (vr) 
meaning respectively, in terni 2 ^! of t defined in (lE.lj) or (|E.2j) 

t = e in(l+iV3)/2 = ie -*V3/2 i (E 3) 

and 

t = e iw(1+i) = -e _7r (E.4) 

that is t 2 = —e~ nv ^ and t 2 = —e~ 2lT respectively. Actually the variable t in (lE.lj) or in [103] is 
exactly what is called the multiplicative crossing in conformal theory [40 [ 1521 [731 180] . Conformal 
field theoreticians are keen on introducing modular group structure for which the multiplicative 
crossing is seen as a modular nome q. If we follow this line recalling the relation q = exp(i7rr) 
between the nome and the half period ratio r (see [100] ). we find that the two previous situations 
actually correspond to singled-out values of the modular j-function namely j((l+iy/3)/2) = (0) 3 



26 Note that this Q, corresponding to the number of state of the Potts model, should not be confused with 
a nome q. It was unfortunately denoted q in |46| . 

27 One has to be careful with the various notations in the literature where, as far as nomes are concerned, one 
moves from q to q 2 . In [100] the nome q corresponds to t 2 . Relation (|9.12|l reads q = q 2 = e 2l7rT = — e - "^. 
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for (|E.3p . and j(l + i) = (12) 3 for (|E,4p . which actually correspond to Heegner numbers and 
complex multiplication [TTJ I1UU] . 

Considering A-extensions of two-point diagonal correlation function of the Ising model, we 
found [25] modular curves corresponding to polynomial relations between a (modular) function 
and its first derivative, this (modular) function being a very simple ratio of Jacobi theta functions 
(see Section 6.1 in [25]). Along this line it is worth recalling the "special value" — e~ nv ^ of the 
nome of Jacobi theta functions (at zero argument) for which a ratio of Jacobi theta functions 
becomes a simple algebraic expression [120] 

<W0,- e -^)_ (W 3_ 7) V*. 



6> 3 (0, -e-^v^) 

At this step it is fundamental to raise an important confusion that overwhelms the theoretical 
physics literature. In many domains of theoretical physics the existence of a modular group 
and/or A-th root of unity situations in some "nome" al way s denoted q, is underlined and 
analyzed. In Liouville theory this nome q is the exponential 2 ^ of H, in conformal field theonl 2 ^! 
two q's, and two modular group structures, can be introduced, the second one corresponding 
to finite size analysis with the introduction of a modular parameter for partition function on 
a (finite size I x /') torus (see for instance (3.33) in [95]). 

Sticking with Baxter's notations the complex multiplication situation, we see in this paper 
with selected values like 1 + 3ui + Aw 2 = 0, corresponds to selected values of the modulus of 
the elliptic curves, or of the nome q which measures the distance to criticality (temperature-like 
variable) of the off-critical lattice model. In contrast the selected values (IB.lj) of A (for which 
modular curves are seen to occur for the A-extensions of the correlation functions) correspond to 
A-th root of unity situations for the multiplicative crossing x. Most of the field theory papers 
(QFT, CFT, . . . ) where selected values (A-th root of unity situations) occur correspond to 
models at criticality: for these models there is no (temperature-like off-critical) variable like 
our previous nome q (the elliptic curve is gone, being replaced by a rational curve). All the 
selected situations encountered are in the multiplicative crossing variable x within a rational 
parametrization of the model. 



F Factorisations of multiple integrals linked to C(3) 

From the series expansion of the triple integral (|9.13p we have obtained the corresponding order 
four Fuchsian linear differential equation (Dx denotes d/dx) 

t n 4 , 2(3x -1) 3 (7x 2 + (n 2 + n - 5)x - 2n(n + 1)) 2 
[X — 1)X [x — l) z x z 

(x 2 + 2n(n + 1)) n(n + 1) ((n 2 + n + l)x + (n - l)(n + 2)) 

{x — l) 2 x 3 (x — 1) 2 X 4 

which has the following factorization in order-one differential operator: 

Ln _ ( Dx + *>sm ( Dx + d J?m ( Dx + ( Dx + 



dx J \ dx J \ dx J \ dx 



28 Not to be confused with the q of the (/-state Potts model in the paper that cope with Liouville theory and 
Potts model in the same time! 

29 They are, of course, many other occurrences of modular groups and/or occurrences of a nome q (quantum 
dilogarithms, g-deformation theories, g-difference equations, g-Painleve, g-analogues of hypergeometric functions, 
. . . ). The confusion is increased with the dilute Al models and their relations with the Ising Model in a Field for 
which the corresponding nome q could be associated with the magnetic field of the Ising model fTHT 11071 11181 1119] . 



38 



S. Boukraa, S. Hassani, J.-M. Maillard and N. Zenine 



where the order-one differential operators have rational solutions since: 

Ay = -{n- l)ln(x) +21n(s- 1) +ln(P„), 
A 2 = (n + 1) ln(x) - (n - 1) ln(x - 1) - ln(P n ) + ln(Q n ), 
A 3 = -nln(x) + (n + 1) ln(x - 1) + ln(P n ) - ln(Q n ), 
= nln(x) — ln(P n ), 

and where P n and Q n are (normalized) polynomials in x of degree n, which satisfy, together 
with Pn and Qn^ (m = 1, . . . , 4), their m-th derivative with respect to x, a system of coupled 
differential equations [30] , 

Such factorization in order-one differential operator having rational solutions is characteristic 
of the strong geometrical interpretation we are seeking for (interpretation of n-fold integrals 
as periods of some algebraic variety) for the Puchsian linear differential operators we have 
obtained for many ?i-fold integrals (of the "Ising class" [8]). Such a factorization in order-one 
linear differential operator having rational solutions does not seem to take place in general for 
our Fuchsian linear differential operators, but seems actually to occur modulo many primes for 
the Fuchsian linear differential operators of the x ■ Such calculations, mixing geometrical 
interpretation and "modular" calculations on our n-fold integrals, remain to be done. 

Acknowledgements 

We have deserved great benefit from discussions on various aspects of this work with F. Chyzak, 
G. Delfino, S. Fischler, P. Flajolet, A.J. Guttmann, M. Harris, I. Jensen, L. Merel, G. Mussardo, 
B. Nickel, J.H.H. Perk, B. Salvy, C.A. Tracy and N. Witte. We thank A. Bostan for a search 
of linear ODEs modulo primes with one of his magma program. We thank one of the three 
referees for very usefull comments. We acknowledge a CNRS/PICS financial support. One of 
us (NZ) would like to acknowledge kind hospitality at the LPTMC where part of this work has 
been completed. One of us (JMM) thanks the MASCOS (Melbourne) where part of this work 
was performed. 

References 

[1] Adler A., Ramanan S., Moduli of Abelian varieties, Lecture Notes in Mathematics, Vol. 1644, Springer- 
Verlag, Heidelberg Berlin, 1996. 

[2] Au-Yang H., Perk J.H.H. , Correlation functions and susceptibility in the Z- invariant Ising model, in Math- 
Phys Odyssee 2001: Integrable Models and Beyond, Editors T. Miwa and M. Kashiwara, Birkhauser, Boston, 
2002, 23-48. 

[3] Au-Yang H., Perk J.H.H., Critical correlations in a Z-invariant inhomogeneous Ising model, Phys. A 144 
(1987), 44-104. 

[4] Au-Yang H., Perk J.H.H., Ising correlations at the critical temperature, Phys. Lett. A 104 (1984), 131-134, 
see equation (4) on page 3. 

[5] Au-Yang H., Perk J.H.H., Star-triangle equations and identities in hypergeometric series, Internat. J. Modern 
Phys. B 16 (2002), 1853-1865. 

[6] Au-Yang H., Perk J.H.H., Wavevector-dependent susceptibility in aperiodic planar Ising models, in Math- 
Phys Odyssee 2001: Integrable Models and Beyond, Editors T. Miwa and M. Kashiwara, Birkhauser, Boston, 
2002, 1-21. 

[7] Babelon O., Bernard D., From form factors to correlation functions: the Ising model, Phys. Lett. B 288 
(1992), 113-120. 

[8] Bailey D.H., Borwein J.M., Crandall R.E., Integrals of the Ising class, J. Phys. A: Math. Gen. 39 (2006), 
12271-12302. 



From Ising Form Factors to n-Fold Integrals 



39 



[9] Ballico E., Casnati G., Fontanari C, On the birational geometry of Moduli spaces of pointed curves, 



math.AG/0701475| 



[10] Barad G., The fundamental group of the real moduli spaces Mo, n - A preliminary report on the topological 
aspects of some real algebraic varieties, http://www.geocities.com/gbarad2002/group.pdf 

[11] Baxter R.J., Exactly solved models in statistical mechanics, Academic Press, London, 1982. 

[12] Baxter R.J., Solvable eight vertex model on an arbitrary planar lattice, Phil. Trans. R. Soc. London A 289 
(1978), 315-346. 

[13] Baxter R.J., Partition function of the eight-vertex lattice model, Ann. Physics 70 (1972), 193-228. 
Baxter R.J., One-dimensional anisotropic Heisenberg chain, Ann. Physics 70 (1972), 323-337. 

[14] Baxter R.J., Eight- vertex model in lattice statistics and one-dimensional anisotropic Heisenberg chain. 

I. Some fundamental eigenvectors, Ann. Physics 76 (1973), 1-24. 

Baxter R.J., Eight- vertex model in lattice statistics and one-dimensional anisotropic Heisenberg chain. 

II. Equivalence to a generalized ice-type lattice model, Ann. Physics 76 (1973), 25-47. 

Baxter R.J., Eight-vertex model in lattice statistics and one-dimensional anisotropic Heisenberg 
chain. HI. Eigenvectors of the transfer matrix and the Hamiltonian, Ann. Physics 76 (1973), 48-71. 

[15] Bazhanov V.V., Mangazeev V.V., The eight-vertex model and Painleve VI, in Special Issue on Painleve VI, 
J. Phys. A: Math. Gen. 39 (2006), 12235-12243, |hep-th/0602122| 

[16] Bazhanov V.V., Nienhuis B., Warnaar O., Lattice Ising model in a field: E$ scattering theory, Phys. Lett. B 
322 (1994), 198-206, [hep-th/9312169l 

[17] Bertin J., Peters C, Variations de structure de Hodge, Varietes de Calabi-Yau et symetrie miroir, Panorama 
et Syntheses, Vol. 3, Societe Mathematique de France, 1996. 

[18] Beukers F., A note on the irrationality of £(2) and £(3), Bull. London Math. Soc. 11 (1979), 268-272. 

[19] Bloch S., Motives associated to graphs, Jpn. J. Math. 2 (2007), 165-196, 

available at http : //math. bu. edu/people/kayeat s/motives/graph_rept 06 1017 .pdf 

[20] Boel R.J., Kasteleyn P.W., Correlation-function identities for general Ising models, Phys. A 93 (1978), 
503-516. 

[21] Boel R.J., Kasteleyn P.W., Correlation-function identities and inequalities for Ising models with pair inter- 
actions, Comm. Math. Phys. 161 (1978), 191-208. 

[22] Boos H.E., Korepin V.E., Evaluation of integrals representing correlations in the XXX Heisenberg spin 
chain, in MathPhys Odyssey 2001, Editors M. Kashiwara and T. Miwa, Birkhfiuser, 2002, 65-108. 

[23] Boos H.E., Gohmann F., Klumper A., Suzuki J., Factorization of multiple integrals representing the density 
matrix of a finite segment of the Heisenberg spin chain, in The 75th Anniversary of the Bethe Ansatz, 
J. Stat. Mech. Theory Exp. 2006 (2006), P04001, 13 pages. 

[24] Boos H.E., Korepin V.E., Quantum spin chains and Riemann zeta function with odd arguments, J. Phys. A: 
Math. Gen. 34 (2001), 5311-5316, [hep-th/0104008| 

[25] Boukraa S., Hassani S., Maillard J.-M., McCoy B.M., Orrick W.P., Zenine N., Holonomy of the Ising model 
form factors, J. Phys. A: Math. Theor. 40 (2007), 75-111, |math-ph/06090"74| 

[26] Boukraa S., Hassani S., Maillard J.-M., McC oy B.M., Weil J.A., Zenine N., Fuchs versus Painleve, J. Phys. A: 
Math. Theor. 40 (2007), 12589-12605, |math-ph/0701014| 

[27] Boukraa S., Hassani S., Maillard J.-M., McCoy B.M., Weil J.A., Zenine N., Painleve versus Fuchs, J. Phys. A: 
Math. Gen. 39 (2006), 12245-12263, |math-ph/0602010| 

[28] Boukraa S., Hassani S., Maillard J.-M., McCoy B.M., Zenine N., The diagonal Ising susceptibility, J. Phys. A: 
Math. Theor. 40 (2007), 8219-8236, |math-ph/0703009| 

[29] Boukraa S., Hassani S., Maillard J.M., Zenine N., Landau singularities and singularities of holonomic inte- 
grals of the Ising class, J. Phys. A: Math. Theor. 40 (2007), 2583-2614, |math-ph/0701016| 

[30] Boukraa S., Hassani S., Maillard J.M., Zenine N., Singularities of n-fold integrals of the Ising class and the 
theory of elliptic curves, J. Phys. A: Math. Theor. 40 (2007), 11713-11748, larXiv:0706.3367l 

[31] Brown F.C.S., Periodes des espaces des modules Mo, n et valeurs zetas multiples. Multiple zeta values and 
periods of moduli spaces M ,„, CPAS C. R. Acad. Sci. Paris Ser I 342 (2006), 949-954. 

[32] Buff X., Fehrenbach J., Lochak P., Schnepps L., Vogel P., Espaces de modules des courbes, groupes modu- 
laires et theorie des champs, Panorama et Syntheses, Numero 7, Societe Mathematique de France, 1999. 

[33] Cecotti S., Vafa C, Ising model and n = 2 supersymmetric theories, Comm. Math. Phys. 157 (1993), 
139-178, |hep-th/9209085| 



40 



S. Boukraa, S. Hassani, J.-M. Maillard and N. Zenine 



Cresson J., Fischler S., Rivoal T., Series hypergeometriques multiples et polyzetas, |math.NT /0609743 

Eden R.J., Landshoff P.V., Olive D.I., Polkinghorne J.C., The analytic S-matrix, Cambridge University 
Press, 1966. 

Erdeleyi A., Asymptotic series, Dover Publishing Co., New York, 1956, p. 47. 

Erdeleyi, Bateman manuscript project, higher transcendental functions, McGraw Hill, New York, 1955. 
Farkas C, Guibney A., The Mori cones of moduli spaces of pointed curves of small genus, Trans. Amer. 



Math. Soc. 355 (2003), 1183-1199, [math.AG/0111268 



Felder G., Varchen ko A., The elliptic Gamma function and SL(3,Z) x Z 3 , Adv. Math. 156 (2000), 44-76, 
|math.QA/9907061| 

Feverati G., Grinza P., Integrals of motion from TBA and lattice-conformal dictionary, Nuclear Phys. B 
702 (2004), 495-515, [hep-th/0405110l 

Fischler S., Groupes de Rhin- Viola et integrales multiples, J. Theor. Nombres Bordeaux 15 (2003), 479-534. 



Fischler S., Integrales de Brown et de Rhin- Viola pour £(3), math.NT/0609799 



Fischler S., Irrationalite de valeurs de zeta (d'apres Apery, Rivoal, . . . ), Seminaire Bourbaki. Expose, no. 910, 
2003, Asterisque 294 (2004), 27-62, |math.NT/0303066| 

Gamier R., Sur les singularites irregulieres des equations differentielles lineaires, J. Math. Pures et Appl. 2 
(1919), 99-198. 

Gerkmann R., Relative rigid cohomology and deformation of hypersurfaces, Int. Math. Res. Pap. IMRP 
2007 (2007), no. 1, Art. ID rpm003, 67 pages. 

Gervais J.-L., Neveu A., Non-standard 2D critical statistical models from Liouville theory, Nuclear Phys. B 
257 (1985), 59-76. 

Glutsuk A. A., Stokes operators via limit monodromy of generic perturbation, J. Dynam. Control Systems 
5 (1999), 101-135. 

Goncharov A.B., Manin Y., Multiple ("-motives and moduli spaces Mo, n , Comp. Math. 140 (2004), 1-14, 
math.AG/0204102| 



Groeneveld J., Boel R.J., Kasteleyn P.W., Correlation function identities for general planar Ising systems, 
Phys. A 93 (1978), 138-154. 

Guzzetti D., The elliptic representation of the general Painleve VI equation, Comm. Pure Appl. Math. 55 



(2002), 1280-1363, math.CV/0108073 



Hanna M., The modular equations, Proc. London Math. Soc. 28 (1928), 46-52. 

Hara Y., Jimbo M., Konno H., Odake S., Shiraishi J., Free field approach to the dilute AL models, J. Math. 
Phys. 40 (1999), 3791-3826. 

Harris M., Potential automorphy of odd-dimensional symmetric powers of elliptic curves, and applications, 
in Algebra, Arithmetic and Geometry - Manin Festschrift, Progress in Mathematics, Birkhauser, to appear. 

Hassett B., Tschinkel Y., On the effective cone of the moduli space of pointed rational curves, in Topology 
and Geometry: Commemorating SISTAG, Contemp. Math. 314 (2002), 83-96, |math.AG/0110231~| 



Huttner M., Constructible sets of linear differential equations and effective rational approximations of 
polylogarithmic functions, Israel J. Math. 153 (2006), 1-44. 

Huttner M., Equations differientielles fuchsiennes; Approximations du dilogarithme, de £(2) et £(3), Publ. 
IRMA, Lille, 1997. 

Jaekel M.T., Maillard J.-M., Inverse functional relations and disorder solutions on the Potts models, 
J. Phys. A: Math. Gen. 17 (1984), 2079-2094. 

Jaekel M.T., Maillard J.-M., Inverse functional relations on the Potts model, J. Phys. A: Math. Gen. 15 
(1982), 2241-2257. 

Jimbo M., Kedem R., Konno H., Miwa T., Weston R., Difference equations in spin chains with a boundary, 
Nuclear Phys. B 448 (1995), 429-456, |he"p-th/9502060| 

Jimbo M., Miwa T., Studies on holonomic quantum fields. XVII, Proc. Japan Acad. Ser. A Math. Sci. 56 
(1980), 405-410, Erratum, Proc. Japan Acad. Ser. A Math. Sci. 57 (1981), 347. 

Jimbo M., Miwa T., Nakayashiki A., Difference equations for the correlations of the eight vertex model, 
J. Phys. A: Math. Gen. 26 (1993), 2199-2209, [hep-th/9211066| 



From Ising Form Factors to n-Fold Integrals 



41 



Katz N.M., Introduction aux travaux recents de Dwork, in Proc. Sympos. Pure Math., Vol. XX (State Univ. 
New York, Stony Brook, New York, 1969), Amer. Math. Soc, Providence, R.I. 1971, 65-75. 

Katz N.M., Nilpotent connections and the monodromy theorem: applications of a result of Turrittin, Pub- 
lications mathematiques de VIHES 39 (1970), 175-232. 

Katz N.M., On the differential equations satisfied by period matrices, Inst. Hautes Etudes Sci. Publ. Math. 
35 (1968), 223-258. 

Katz N.M., Rigid local systems, Ann. of Math. Stud., Vol. 139, Princeton University press, 1996. 

Katz N.M., Travaux de Dwork, in Seminaire Bourbaki (1971/1972), Exp. No. 409, Lecture Notes in Math., 
Vol. 317, Springer Verlag, 1973, 167-200. 

Kaufman B., Crystal statistics. II. Partition function evaluated by spinor analysis, Phys. Rev. 76 (1949), 
1232-1243. 

Kaufman B., Onsager L., Short-range order in a binary Ising lattice, Phys. Rev. 76 (1949), 1244-1252. 

Krattenthaler C, Rivoal T., An identity of Andrews, multiple integrals, and very- well-poised hypergeometric 
series, Ramanujan J. 13 (2007), 203-219, [m"ath.CA/0312l48| 

Kreimer D., Knots and Feynman diagrams, Cambridge Lecture Notes in Physics, Vol. 13, Cambridge Uni- 
versity Press, 2000, Chapter 9. 

Lario J.-C, Elliptic curves with CM defined over extensions of type (2,..., 2), available at 
http : //www-ma2 .upc . es/~lario/ellipticm.htm 

Levin A., Racinet C, Towards multiple elliptic polylogarithms, |math. NT /0703237 

Lukyanov S., Pugai Y., Multi-point local height probabilities in the integrable RSOS model, Nuclear Phys. B 
473 (1996), 631-658, |hep-th/9602074l 

Lyberg I., McCoy B.M., Form factor expansion of the row and diagonal correlation functions of the two 
dimensional using model, J. Phys. A: Math. Theor. 40 (2007), 3329-3346, |math-ph/0612051| 

Maier R.S., On rationally parametrized modular equations, math. NT/0611041 



Maillard J.-M., The inversion relation: some simple examples, J. Physique 46 (1984), 329-341. 

Maillard J.-M., Boukraa S., Modular invariance in lattice statistical mechanics, Ann. Fond. Louis de Broglie 
26 (2001), Special Issue 2, 287-328. 

Manin Yu.L, Sixth Painleve equation, Universal elliptic curve, and mirror of P 2 , Amer. Math. Soc. Transl. 
Ser. 2 186 (1998), 131-151, [alg-geom/96050lo| 



Manojlovic N., Nagy Z., Creation operators and algebraic Bethe ansatz for elliptic quantum group E TiV (so3) 



J. Phys. A: Math. Theor. 40 (2007), 4181-4191, math.QA/0612087 



Martinez J.R., Correlation functions for the ^-invariant Ising model, hep-th/9609135 



Mazzocco M., Picard and Chazy solutions to the Painleve VI equation, Math. Ann. 321 (2001), 157-195, 
|math.AG/9901054| 

McCoy B.M., Perk J.H.H., Wu T.T., Ising field theory: quadratic difference equations for the n-point 
Green's functions on the lattice, Phys. Rev. Lett. 46 (1981), 757-760. 

McCoy B., Tracy OA., Wu T.T., Painleve equations of the third kind, J. Math. Phys. 18 (1977), 1058-1092. 

McCoy B.M., Wu T.T., Nonlinear partial difference equations for the two-dimensional Ising model, Phys. 
Rev. Lett. 45 (1980), 675-678. 

Murata M., Sakai H., Yoneda J., Riccati solutions of discrete Painleve equations with Weyl group symmetry 
of type J. Math. Phys. 44 (2003), 1396-1414, rnlin.SI/0210040| 

Mussardo G., II Modello di Ising introduzione alia teoria dei campi e delle transizioni di fase, Editor Bollati 
Boringhieri, 2007. 

Nickel B., On the singularity structure of the Ising model susceptibility, J. Phys. A: Math. Gen. 32 (1999), 
3889-3906. 

Nickel B., Addendum to 'On the singularity structure of the Ising model susceptibility', J. Phys. A: Math. 
Gen. 33 (2000), 1693-1711. 

Nickel B., Comment on "The Fuchsian differential equation of the square lattice Ising model x' 3 ^ susceptibi- 
lity" [J. Phys. A: Math. Gen. 37 (2004), 9651-9668] by N. Zenine, S. Boukraa, S. Hassani and J.-M. Maillard, 
J. Phys. A: Math. Gen. 38 (2005), 4517-4518. 



42 



S. Boukraa, S. Hassani, J.-M. Maillard and N. Zenine 



[90] Okamoto K., Studies on the Painleve equations. I. Sixth Painleve equation P VI, Ann. Mat. Pura Appl. (4) 
146 (1987), 337-381. 

[91] Onsager L., Crystal statistics. I. A two-dimensional model with an order disorder transition, Phys. Rev. 65 
(1944), 117-149. 

[92] Onsager L., Nuovo G'imento 6 (1949), suppl., 261. 

[93] Orrick W.P., Nickel B.G., Guttmann A. J., Perk J.H.H., The susceptibility ol the square lattice Ising model: 
new developments, J. Statist. Phys. 102 (2001), 795-841, [cond-mat/0103074| 

[94] Palmer J., Tracy C.A., Two-dimensional Ising correlations: convergence of the scaling limit, Adv. in Appl. 
Math. 2 (1981), 329-388. 

[95] Pearce P.A., Temperley-Lieb operators and critical A-D-E models, Internat. J. Modern Phys. A 4 (1990), 
715-734. 

[96] Perk J.H.H., Quadratic identities for Ising model correlations, Phys. Lett. A 79 (1980), 3-5. 

[97] Perk J.H.H., Au-Yang H., Some recent results on pair correlation functions and susceptibilities in exactly 
solvable models, in Dunk Island Conference in Honor of 60th Birthday of A.J. Guttmann, J. Phys. Con]. 
Ser. 42 (2006), 231-238, [math-ph/0606046l 

[98] Perk J.H.H., Capel H.W., Time-dependent ia-correlation functions in the one-dimensional A"Y-model, 
Phys. A 89 (1977), 265-303, see equation (6.16). 

[99] Picard E., Memoire sur la theorie des functions algebriques de deux variables, Journal de Liouville 5 (1889), 
135-319. 



[100 
[101 

[102; 

[103 
[104 
[105; 

[106 
[107 

[108; 

[109 
[110 
[111 
[112 
[113 
[114 
[115; 



Piezas T. Ill, Weisstein E.W., j-function, from MathWorld-A Wolfram Web Resource, 
http : //mathworld. wolf ram. com/ j -Funct ion . html | 

Racinet G., Doubles melanges des polylogarithmes multiples aux racines de l'unite, Publ. Math. Inst. Hautes 
Etudes Sci. No. 95 (2002), 185-231, [math.QA/0202142| 

Ramis J.-P., Confluence et Resurgence, J. Fac. Set. Univ. Tokyo Sect. IA Math. 36 (1992), 703-716. 

Rammal R., Maillard J.-M., g-state Potts model on the checkerboard lattice, J. Phys. A: Math. Gen. 16 
(1983), 1073-1081. 

Rammal R., Maillard J.-M., Some analytical consequences of the inverse relation for the Potts model, 
J. Phys. A: Math. Gen. 16 (1983), 353-367. 

Rhoades R.C., Elliptic curves and modular forms (notes based on A Course at the University of Wis- 
consin - Madison MATH 844 during the Spring 2006 taught by Professor Nigel Boston), available at 
http : //www. math, wise . edu/~rhoades/Notes/EC .pdf 

Sakai H., Rational surfaces associated with affine root systems and geometry of the Painleve equations, 
Comm. Math. Phys. 220 (2001), 165-229. 

Seaton K.A., Batchelor M.T., The dilute A4 models, the £7 mass spectrum and the tricritical Ising model, 
J. Math. Phys. 43 (2002), 2636-2653, |math-ph/0110021| 

Singer M.F., Testing reducibility of linear differential operators: a group theoretic perspective, Appl. Alg. 
Eng. Commun. Comp. 7 (1996), no. 2, 77-104. 

Sorokin V.N., On the measure of transcendency of the number 7r 2 , Mat. Sb. 187 (1996), no. 12, 87-120 
(English transl: 56. Math. 187 (1996), 1819-1852). 

Sorokin V.N., Apery's theorem, Vestnik Moskov. Univ. Ser. I Mat. Mekh. (1998), no. 3, 48-53, 74 (English 
transl.: Moscow Univ. Math. Bull. 53 (1998), no. 3, 48-52). 

Stanev Y.S., Todorov I., On the Schwartz problem for the SU2 Knizhnik-Zamolodchikov equation, Lett. 
Math. Phys. 35 (1995), 123-134. 

Sternin B.Yu., Shatalov V.E., On the confluence phenomenon of fuchsian equations, J. Dynam. Control 
Systems 3 (1997), 433-448. 

Todorov I.T., Arithmetic features of rational conformal field theory, Ann. Inst. H. Poincare 63 (1995), 
427-453. 

van der Put M., Singer M.F., Galois theory of linear differential equations, Grundlehren der Mathematischen 
Wissenschaften, Vol. 328, Springer- Verlag, Berlin, 2003, available at http://www4.ncsu.edu/singer/ 

van Hoeij M., Rational solutions of the mixed differential equation and its application to factorization of 
differential operators, in Proceedings ISSAC '96, ACM, New York, 1996, 219-225. 



From Ising Form Factors to n-Fold Integrals 



43 



Vasilyev D.V., On small linear forms for the values of the Riemann zeta-function at off integers, Doklady 
NAN Belarusi (Reports of the Belarus National Academy of Sciences) 45 (2001), no. 5, 36-40 (in Russian). 

Vasilyev D.V., Some formulas for Riemann zeta-function at integer points, Vestnik Moskov. Univ. Ser. I 
Mat. Mekh. (1996), no. 1, 81-84 (English transl: Moscow Univ. Math. Bull. 51 (1996), no. 1, 41-43). 

Warnaar O., Nienhuis B., Seaton K.A., New construction of solvable lattice models including an Ising model 
in a field, Phys. Rev. Lett. 69 (1992), 710-712. 

Warnaar O., Pearce P.A., Exceptional structure of the dilute A3 model: Eg and E7 Rogers-Ramanujan 
identities, J. Phys. A: Math. Gen. 27 (1994), L891-L897, |hep-th/9408136| 

Weisstein E.W., Jacobi theta functions, from MathWorld-A Wolfram Web Resource, 
http : //mathworld. wolf ram. com/ JacobiThetaFunctions .html 

Witte N.S., Isomonodromic deformation theory and the next-to-diagonal correlations of the anisotropic 
square lattice Ising model, J. Phys. A: Math. Theor. 40 (2007), F491-F501. l arXiv:0705.0557l 

Wu T.T., Theory of Toeplitz determinants and of the spin correlations of the two-dimensional Ising model, 
Phys. Rev. 149 (1966), 380-440. 

Wu T.T., McCoy B.M., Tracy C.A., Barouch E., Spin-spin correlation functions for the two dimensional 
Ising model: exact theory in the scaling region, Phys. Rev. B 13 (1976), 316-374. 

Yamada K., On the spin-spin correlation function of the Ising square lattice and the zero field susceptibility, 
Progr. Theoret. Phys. 71 (1984), 1416-1418. 

Yang C.N., The spontaneous magnetization of the two dimensional Ising model, Physical Rev. (2) 85 (1952), 
808-816. 

Zenine N., Boukraa S., Hassani S., Maillard J.M., The Fuchsian differential equation of the square Ising 
model x (3) susceptibility, J. Phys. A: Math. Gen. 37 (2004), 9651-9668, |math-ph/0407060| 

Zenine N., Boukraa S., Hassani S., Maillard J.M., Ising model susceptibility: Fuchsian differential equation 
for x (4) and its factorization properties, J. Phys. A: Math. Gen. 38 (2005), 4149-4173, |cond-mat/0502155| 

Zenine N., Boukraa S., Hassani S., Maillard J.M., Square lattice Ising model susceptibility: connection ma- 
trices and singular behavior of x (3> and x (4> , J- Phys. A: Math. Gen. 38 (2005), 9439-9474, |hep-th/0506214[ 



math-ph/0506065 



Zenine N., Boukraa S., Hassani S., Maillard J.M., Square lattice Ising model susceptibility: series expansion 
method, and differential equation for x <3) , J- Phys. A: Math. Gen. 38 (2005), 1875-1899, |hep-ph/04lT05l"j 

Zhang O, Confluence et phenomenes de Stokes, J. Math. Sci. Univ. Tokyo 3 (1996), 91-107. 



